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A new metastable isomer of gold was found. The energy 
level system of gold has been studied in the region from the 
threshold to 3.3 Mev. The energy of the metastable state 
was found to be approximately 250 kev and higher activa- 
tion levels which combine with the metastable state are 
found at 1.22+0.03, 1.68, 2.15, 2.56, and 2.97 Mev, re- 
spectively. This excitation has also been produced by fast 


neutrons by the Au-n-n reaction and the threshold for this 
process has been determined. The previously known 43- 
minute period of mercury has been produced by direct 
x-ray irradiation as well as a 42-day period in columbium. 
The strontium activity, known previously from neutron 
excitation, of 2.7 hours half-life has been produced by both 


x-ray and electron bombardment of this element. 





INTRODUCTION 


LTHOUGH activity may be produced in 
stable nuclei by means of the ordinary re- 
actions such as y-m, n-2n, n-y, n-a, p-n, etc., the 
most direct proof of the existence of a metastable 
state in such stable nuclei lies in the production 
of this activity by the direct irradiation of the 
element by x-rays or electrons below the thresh- 
old for photo-disintegration of the given nucleus. 
At these energies the only activity which can be 
produced is that which is excited in one of the 
isotopes originally present and, therefore, can be 
attributed to a metastable state of a stable 
nucleus. If the metastable state belongs to a 
relatively abundant isotope, one can produce 
sufficient activity, even at low energies, and thus 
obtain an excitation curve! from which may be 
obtained much valuable information concerning 
the nuclear energy level system of that nucleus.’ 
The Van de Graaff-Herb electrostatic generator 


1B. Waldman, G. B. Collins, and E. M. Stubblefeld, 
Phys. Rev. A55, 1129 (1939); M. L. Wiedenbeck, Phys. 
Rev. 67, 92 (1945). 
*E. Guth, Phys. Rev. 59, 325 (1941). 





furnishes a convenient and suitable source of 
x-rays and electrons for the study of these 
excitation curves. In a quantitative determina- 
tion of the energies of the various activation 
levels which combine with the metastable state, 
the x-ray excitation curve is the more accurate 
of the two. When the energy of the monokinetic 
electron beam, which produces the x-rays, reaches 
a higher activation level, the curve shows a sharp 
increase in slope. The intersection of the two 
linear portions of the curve, each obtained from 
numerous experimental points, clearly defines the 
energy of this activation level. It is not possible to 
produce excitation by a natural source of y-rays 
since the process is one of line absorption, the 
natural width of the line being of the order of 
millivolts while the total width, including the 
Doppler effect, is of the order of a few volts. 


EXPERIMENTAL 
A. Mercury 


An activity with a half-life of 43 minutes has 
been observed and chemically identified as mer- 
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Fic. 1. Thin walled aluminum counter. The bombarded 
mercury was contained between the outer cylinder and the 
counter wall. 


cury by Heyn,* and others*® when mercury is 
irradiated with fast neutrons in the process 
Hg-n-2n. This activity was also produced by 
Krishnan and Nahum‘ by the Hg-d- reaction. It 
was not clear whether this activity was to be 
ascribed to a metastable state of a stable mercury 
nucleus or to an unstable nucleus which decays to 
Au by K-capture.® 

We have now produced the activity with this 
period by bombarding mercury with x-rays of 
maximum energy of 2.6 Mev. The metallic mer- 
cury was irradiated for thirty minutes with the 
x-rays produced when an electron beam current 
of 150 microamperes struck a thick gold target. 
The mercury was then poured around a thin 
walled aluminum counter (Fig. 1). The decay 
curve shown in Fig. 2 clearly indicates the 
presence of the 43-minute period which can, 
therefore, definitely be attributed to Hg*. The 
relatively small activity prohibits the obtaining 
of a complete excitation curve at the present 
time and indicates that the isotope possessing the 
metastable state is relatively low in abundance. 


B. Strontium 


The isotope Sr*’ has been shown to possess a 
metastable level 0.38 Mev above the ground 
state. This has been excited by the Sr*’-n-n%® 
reaction as well as having been grown in the 
processes: Rb-p-n,° Sr-n-y,** Zr-n-a,!® and from 


*F,. Heyn, Nature 139, 842 (1937). 
4E. McMillan, M. Kamen, and S. Ruben, Phys. Rev. 52, 
375 (1932). 
5M. L. Pool, J. M. Cork, and R. L. Thornton, Phys. Rev. 
52, 239 (1937). 
*R. S. Krishnan and E. A. Nahum, Proc. Camb. Phil. 
Soc. 36, 490 (1940). 
046) A. Dubridge and J. Marshall, Phys. Rev. 57, 348 
1940). 
8 H. Reddemann, Naturwiss. 28, 110 (1940). 
(1935) A. Dubridge and J. Marshall, Phys. Rev. 56, 706 
10 R. Sagane, S. Kojima, G. Miyamoto, and M. Ikawa, 
Phys. Rev. 57, 1180 (1940). 
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Y*’ by K decay.’ This 2.7-hour activity has now 
also been produced by the direct x-ray bombard- 
ment of strontium in the reaction : Sr*?+y—Sr87* 
and by the electron bombardment of strontium 
chloride in the reaction: 


Sr®7+¢e,;—Sr®7* +e,>. 


C. Columbium 


A radioactive isomer of stable Cb® has been 
produced in the beta-decay of Zr®.!° Sagane, 
et al.!° have reported a half-life of the order of 55 
days. Pool and Edwards" have more recently 
reported a decay period of 38.7 days for this 
activity. After eighteen hours of bombardment 
with x-rays produced by a beam of 150 micro- 
amperes of 2.6-Mev electrons stopped on a gold 
target, a small columbium sample showed an 
initial activity of 2000 counts per hour when the 
sample was placed before a thin windowed (1.5- 
mil Al) counter. This activity was found to decay 
with a half-life of approximately 42 days. 
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Fic. 2. Decay curve for Hg*. The half-life of mercury in the 
transition: Hg*—+Hg+-7 is 43 minutes. 


1 Q. Hahn and F. Strassmann, Zeits. f. Physik 121, 729 
(1943). 
be W. N. Moquin and M. L. Pool, Phys. Rev. 65, 60 
1944). 

%M. L. Pool and J. E. Edwards, Phys. Rev. 67, 00 
(1945). 
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D. Gold 


We have previously reported" an activity of 
half-life 7.50.5 seconds which was produced by 
the bombardment of gold with an x-ray beam. 
Since there is but one stable isotope of gold, this 
activity must be attributed to a metastable state 
in Au!®?, 

In the ground state gold has a spin of 3/2 
which, therefore, would imply a spin of 9/2 or 
11/2 for the metastable state. Gold is the first 
known isotope in which the spin in the ground 
state is other than 1/2 or 9/2 and is thus an 
exception to the empirical rule set forth by 
Mattauch"® that stable isotopes possessing such 
metastable states will have spins of 1/2 or 29/2. 

This activity has been determined as a function 
of the accelerating potential from the threshold to 
3.3 Mev and is shown in Fig. 3. The activity is 
given as the number of counts during the first 
minute starting five seconds after the sample had 
been bombarded for thirty seconds with x-rays 
produced by an electron beam of 100 micro- 
amperes. The sample which was used in these 
measurements was a gold cylinder which formed 
the cathode of a self-quenching argon-ether 
counter. Several of such cylinders had been tried 
including a copper cylinder plated on the inside 
with $ mil of gold. However, a pure gold cylinder 
was found to be most satisfactory and was used 
in obtaining all measurements. 

As in the previous cases, it is seen that the 
excitation curve is composed of a series of straight 
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Fic. 3. X-ray excitation curve for gold. The intersection 
of the straight line segments gives the energy of the 
activated levels which combine with the metastable state. 
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16 J. Mattauch, Zeits. f. Physik 117, 246 (1941). 


energy 
mev 


2.97 


Spin 





2.56 





2.15 





1.68 





1.22 




















0.25 f y y y 
@) y 


























3/2 


Fic. 4. Nuclear energy level diagram for gold showing the 
various transitions which occur. 


line segments, the intersection of which gives the 
positions of the higher activation levels which 
combine with the metastable state. These breaks 
are found at 1.22+0.03 Mev, 1.68, 2.15, 2.56 and 
2.97 Mev, respectively. 

The energy of the metastable level was ob- 
tained by measuring the energy of the conversion 
electrons emitted in the transition : Au'®7*—>Au!®” 
++. A thick gold disk was placed before the thin 
window of a counter with a brass cylinder, and the 
sample was irradiated at a fixed accelerating 
potential and with a constant beam current. 
Aluminum foils were interposed between the 
counter and the gold, and thus the activity was 
determined as a function of the thickness of 
aluminum. From this data, one obtains a value of 
approximately 250 kev for the energy of this 
level. 

A nuclear energy level diagram showing the 
various levels which combine with the metastable 
level and the transition involved is shown in 
Fig. 4. 

Activities due to metastable states of stable 
nuclei have been produced in several instances? !* 
by the inelastic collision of fast neutrons with the 


16 M. Goldhaber, R. D. Hill, and L. Szilard, Phys. Rev. 
55, 47 (1939); Nature 142, 521 (1938). 
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Fic. 5. Target arrangement for the excitation of gold with 
neutrons produced from the Be-y-n reaction. 


stable nucleus. Up to this time no quantitative 
data have been obtained concerning the energy 
of the neutrons required to produce this excita- 
tion. One should expect, of course, that the 
neutron threshold is practically identical with the 
x-ray threshold, but a direct proof is of interest. 
The electrostatic generator furnishes a convenient 
source of neutrons of known energy which can be 
produced by the Be-y-n reaction. The distribu- 
tion of energies of the neutrons is continuous as is 
the x-ray spectrum which produces them. 

The intensity of any neutron line is linear!’ 
(in this energy region) with the applied ac- 
celerating potential above its threshold. The 
maximum energy of the neutrons so obtained is: 
V,=8/9(V—1.63), where V is the applied elec- 
tron accelerating voltage. 

Since the threshold for the gold excitation was 
found to be 1.22 Mev with x-rays, it seemed 
possible that this activity could also be produced 
with the neutron spectrum described above. 
Gold has the additional advantage over any 
other known metastable isomer since it exhibits 
no other short period neutron capture reactions. 

An arrangement was set up as shown in Fig. 5. 
The x-rays were allowed to strike one-half pound 
of powdered beryllium placed directly in front of 
the target. Lead of four-inch thickness was placed 
between the beryllium and the self-quenching 
counter containing the gold cathode. The lead 
served as a shield to protect the counter from the 
x-rays. 

With the apparatus set up in this manner, the 


17M. L. Wiedenbeck and Rev. C. J. Marhoefer, Phys. 
Rev. 67, 54 (1945). 
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activity produced by the x-rays was determined 
by bombardment without having the beryllium 
in position. Next, the beryllium was placed in 
position. No increase in activity beyond the small 
residual x-ray activity was found until the 
maximum energy of the neutrons exceeded 1.22 
Mev. With maximum neutron energies in the 
range from 1.22 to 1.3 Mev, a considerable excess 
counting rate was obtained which increased 
linearly with energy as in the case of x-ray excita- 
tion. This indicates that the threshold for the 
process, Au!®?-+-2,,;—>Au!®’* +-n,», is identical with 
that of x-rays. 


DISCUSSION 


Again, as in the case of silver, cadmium, and 
indium, the energy levels of gold were found to 
start in the region of 1.2 Mev, and successive 
higher levels are approximately equally spaced 
with a difference of 0.4 Mev between them. No 
convergence of the levels is indicated at these 
energies. 

These experiments check the theory and show 
that the excitation processes in the cases of x-rays, 
electrons, and neutrons are similar. As one in- 
creases the energy of the x-rays or electrons 
producing the excitation, one eventually reaches 
a limit when this energy is equal to the binding 
energy of a neutron in the nucleus and thus 
produces the photo- or electro-disintegration of 
that nucleus. The similarity with the neutron 
process indicates that when the energy of the 
incident neutrons reaches the binding energy of a 
nuclear neutron, the m-2n process can occur.!* 
This process may be looked upon as one of neutro- 
disintegration, differing from that of photo- 
disintegration, however, due to the required 
momentum consideration of the incident neutron. 
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A new theoretical mixture rule is developed for determining the refractive index of colloidal 

. particles from differential refractometric measurements. The new rule gives better results for 
colloidal solutions than the mixture rules of Newton, Lorenz-Lorentz, or Wiener if the difference 
between the refractive indices of the colloidal particles and of the medium is small. If the differ- 
ence is large, but not in excess of 0.7, the refractive index of the particles can still be obtained 
with an accuracy in the third or fourth decimal by using an empirical correction equation. If 
colloidal solutions are very opaque, the refractive index of the particles can be calculated from 
light scattering. Application of the respective equation requires that the refractive index of the 
medium and the size of the particles are known. The gravimetric and volumetric measurements 
involved in the determination of density of colloidal particles can be replaced in part or entirely 


by measurements of light scattering and refractivity. 





I. A NEW MIXTURE RULE FOR DETERMINING THE 
REFRACTIVE INDEX OF COLLOIDAL 
CRYSTALS OR MOLECULES 


A. The Classical Mixture Rules 


SERIES of mixture rules has been pro- 
posed for determining the refractive index 
of a dispersed substance, m,, indirectly, from 
measurements of the refractive index of the total 
dispersed system, ;, and of the refractive index 
of the medium of dispersion, ,,. Except for a few 
rules of minor significance, they are: 
Newton, theoretical : 
n?= PmNm? + pny’, (1) 
Lorenz-Lorentz, theoretical : 


= Om + 
ni?+2 Nm? +2 


n,?—1 


Pp ’ 
n,*+2 








(2) 


Wiener, theoretical : 





Ny? — Nn? Ny? — Ny? 
—memnene GS Ghytrenennasem, (3) 
n?+u Ny? +u 
Lichtenecker, empirical : 
log 2:= om log Im+ ¢> log np, (4) 


Avogadro-Biot-Beer-Landolt-Christiansen- 
Wintgen, empirical : 


N= PmNmt PpMy. (S) 


¢m and g, are the respective volume fractions. 








The term u in Eq. (3) is a function of the shape 
of the particles. u=2n,,, if the particles are 
spherical. Equation (3) is then identical with 
Eq. (2). 

A large series of experiments by various au- 
thors has shown that the “refractive increment” 
(n,—MNm)/c—in which c is the weight concentra- 
tion of the particles—is practically independent 
of ¢ in dilute colloidal solutions. A priori, the 
empirical Eq. (5) appears therefore to be most 
suited for colloidal systems. All equations are 
tested in Table I in which n, is calculated for 
colloidal sulphur." 

In addition to being empirical, Eq. (5) has a 
second disadvantage, which is shared also by the 
other mixture rules: both m, and m, must be 
known. If they are determined in two separate 
experiments, small temperature differences will 
play a role, because they affect the values of 
Nt, Mm, and ¢. 





TABLE I. Refractive indices of S at room temperature. 








Experimen- 
tal values 


Values calculated for infinite dilution from mixture rules 
Assumed density: 2.000) 


Eq. (1) Eqs. (2,3) Eq. (4) Ea. (5) Eq. (9) 





2.085 1.885 2.267 2.184 2.008 1.992 








_} Since sulphur is only slightly lyophilic, the medium of 
dispersion should have little effect upon mz, i.e., it should be 
expected that Np exp =>» cale- 

























It appeared desirable, therefore, to develop a 
new mixture rule which has a theoretical founda- 
tion and which allows the direct use of differential 
refractometric data. 


B. The New Mixture Rule 


A mixture rule which meets the two require- 
ments can be derived easily on the basis of the 
two principal theories developed on light-scat- 
tering. There is, on the one hand, the well-known 
theory of Rayleigh? and Mie,’ in which light 
scattering is considered as an effect due to 
material units. As a consequence, ”, enters into 
the equations developed, in addition to ,,. On the 
other hand, there is the theory of Smoluchowski* 
and Einstein’ in which light scattering is con- 
sidered as an effect due to fluctuations, fluctua- 
tions in density in liquids and fluctuations in 
concentration in liquid mixtures. More complete 
equations for liquid mixtures and solutions have 
been derived by Raman® who takes into account 
in these systems both fluctuations in density and 
in concentration. Light scattering due to fluctua- 
tions in density is negligible compared to light 
scattering due to fluctuations in concentration, if 
the solutions are colloidal. This special case has 
been considered recently by Debye.’ In the equa- 
tion at which Debye arrives, the refractive index 
of the particles obviously does not enter, but ; 
enters in addition to mm. 

A mixture rule can therefore be obtained by 
combining the equation of Rayleigh and the 
equation of Debye. This can be achieved by a 
simple transformation of the former equation. 
The classical formulation of the Rayleigh equa- 


tion is 
24x* sa?—1 
R= (- ) yv*I, (6) 

r4 249 


where R is the Tyndall radiation, Ip the intensity 
of the incident primary beam, a=n,/nm, v the 
number of the light scattering particles per cm’, v 
the volume of these particles, and \ the wave- 








2 Lord Rayleigh, Phil. Mag. 41, 107, 274, 447 (1871) ; 44, 
28 (1897); 47, 375 (1899). 

3G. Mie, Ann. d. Physik 25, 377 (1908); P. Debye, 
Thesis, Munich (1908). 

4M. v. Smoluchowski, Ann. d. Physik 25, 205 (1908). 

5 A, Einstein, Ann. d. Physik 33, 1275 (1910). 

*C. V. Raman and K. R. Ramanathan, Phil. Mag. 45, 
213 (1923); see also R. Gans, Zeits. f. Physik 17, 371 (1923). 
7P. Debye, J. App. Phys. 15, 338 (1944). 
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length in the medium. The quantities R and I 
which are inconvenient for the present purpose 
can be replaced by considering rather the decrease 
in intensity of the primary beam during its 
passage through the light scattering system. 
Since the light intensity can be defined as the 
energy incident upon the unit surface area per 
unit time, the absolute loss in intensity (J»—J) 
has the dimension ergXcm~*Xsec.-!=g Xcm™ 
Xsec.~*. Since the light radiation has the dimen- 
sion, ergXsec.~!, it follows that the Tyndall 
radiation, which is proportional to the volume 
illuminated, has also the dimension gXcm™ 
Xsec.—*. Consequently, the absolute loss in in- 
tensity (Ij —I)=R (whenever the decrease in 
intensity of the primary beam is due exclusively 
to coherent light scattering), and the fractional 
decrease in intensity over unit distance: 


Igp—I 24n* 7a*?—1\? 
ra) et (6a) 
Io At \a?+2 


After substituting the wave-length in vacuum, 
Xo, for X and after replacing the ratio, a, by 
differences of the refractive indices involved, the 
transformed equation : 
242? nn*(n,? — Mm)? 
T= gv, (6b) 
Ao! (ny?+21m?)? 
is obtained, which is suitable for the present 
purpose. Debye’s equation, on the other hand, is: 


32x 1 
Nm? (N,— Nm)?—-. (7) 
3X0! Vv 











T= 


Upon combining (6b) and (7) the mixture rule 


Nm? + 2bNm?(Ny— Mm) 
Ns)? = — “ (8) 
Nm — b(ny— Nm) 





is obtained where 
b=2/3¢. 


Compared to Eq. (5), the measurement of ; is 
replaced by a differential measurement (”;—1m) 
under elimination of temperature errors. The 
advantage of the differential measurement may 
be expected to be particularly pronounced in very 
dilute solutions where (;~—m,,) is very small. 
Such solutions are not likely to give reliable 
results unless differential measurements are 
carried out. 
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TABLE II. Magnitude of error in calculations of refractive 
indices from Eq. (8). 











Errors occur 
Np(8) —m np(s) if mm = 1.333000 in the 

0.00235-0.00775 1.33535-1.34075 6th decimal 
0.00775-0.0245 1.34075-1.3575 5th decimal 
0.0245 —0.082 1.3575 -1.415 4th decimal 
0.082 -0.225 1.415  -1.588 3rd decimal 
0.225 -0.750 1.588 -2.083 2nd decimal 
0.750 -2.18 2.083 -3.51 1st decimal 








C. The Limits of Validity of the New 
Mixture Rule 


The theory of Rayleigh-Mie has been verified 
repeatedly,® i.e., particle sizes determined on the 
basis of the equation developed by Mie, of which 
the Rayleigh Eq. (6) is a special case, were found 
to be in quantitative agreement with the actual 
particle sizes. Similarly, the equation of Debye 
has been verified. There can therefore be no 
doubt that the new mixture rule is valid for 
systems which satisfy the assumptions involved 
in the two theories on light scattering: Eq. (6) is 
based on the assumption that the difference 
(n»—Mm) is very small; Eq. (7) is based on the 
assumption that fluctuations in density are 
negligible compared to fluctuations in concentra- 
tion. Equation (8) should, therefore, be valid, if 
(n»—Mm) is very small and if the dispersed 
particles (crystals or molecules) are large enough 
to scatter considerably more light than the mole- 
cules of the medium. The former condition can be 
fulfilled by choosing the right medium. The latter 
condition is fulfilled in colloidal solutions. It is 
not fulfilled in liquid mixtures nor in solutions of 
small molecules or ions. 

Both Eqs. (6) and (7) presuppose that the 
dispersed systems are dilute and that the light 
scattering particles are small as compared to the 
wave-length used, that they are spherical and 
isotropic. However, identical correction factors 
would have to be introduced into either equation 
in order to make them applicable to systems of 
moderate concentration and to large or non- 
spherical or intrinsically anisotropic particles. 
Consequently, these partly uriknown correction 
factors are cancelled out in the combination of 
Eqs. (6) and (7) to Eq. (8). No influence of the 


* E.g., Niels Pihlblad, Thesis; Uppsala (1918). 
* W. Heller and H. B. Klevens, Phys. Rev. 67, 61 (1945). 


size or shape of colloidal particles upon m,,s) 
should therefore be expected unless the refractive 
index actually varies with the size and shape of 
the particles, an experimentally unknown effect 
which might be found, however, in a spectral 
region near an absorption band.!° Neither should 
the colloid concentration have any effect on m»,s) 
as long as the concentration is small enough so as 
not to affect the refractive indices themselves. 
The probable error involved in calculations of 
ns) and as due to the finite value of (n,—m,,) is 
compiled in Table II. (The error is obtained by a 
procedure discussed below.) The error occurs in 
the fifth or fourth decimal, if (mis) —mm) has a 
value not in excess of 0.025 and 0.082, respec- 
tively. 


D. Another Formulation of the New 
Mixture Rule 


Upon reviewing the present paper, Professor 
Debye drew the author’s attention to another 
simple possibility of deriving the result expressed 
by Eq. (8). The same result can be obtained by 
calculating (Clausius-Mosotti) what the dielectric , 
constant is of a medium in which a sufficiently 
small number of spheres of different dielectric 
constant is dispersed. 

According to Debye, the result expressed by 
the new mixture rule can be formulated as 
follows: 


a=A/(1+ $A), (8a) 
or 
A=a/(1—}a). (8b) 
The quantity 
Ny? — Nn? 
4=——__, 
Nm 


and the quantity 
(1: —Mm)/Mm 


¢ 





a= 


It is assumed, as in the case of Eqs. (6) and (7) 
that the dielectric constant «= n*. Equations (8a) 
and (8b) check, of course, with Eq. (8). 


E. Application of the New Mixture Rule Outside 
of the (n,—n,,) Range of Satisfactory Validity 


It’ is sometimes not possible to select the 
medium in which (n,—m,,) is sufficiently small in 


10 W. Heller, Comptes rendus 205, 971 (1937). 








8 WILFRIED HELLER 


TABLE III. Refractive indices in hydrosols. 








(np(s»—p(H) (p(s) — p(s) 
x<10* x<103 





Dispersed substance ¢ X10? ne Nm Np(8) Np(9) Np(5) 
SiO: 3.075 1.33705 1.3325 1.4837 1.4805 1.4803 3.4 0.2 

SiO: 0.5815 1.33353 1.3325 1.5144 1.5097 1.5099 4.5 —0.2 
a 3.10 1.33946 1.3325 1.5651 1.5571 1.5571 8.0 0.0 
latin 3.10 1.34110 1.3332 1.5988 1.5882 1.5881 10.7 0.1 
Hg-sulfosalicylic acid 2.323 1.34388 1.3330 1.8457 1.8014 1.8015 44.2 —0.1 
AIOOH 1.047 1.33789 1.3325 1.9035 1.8472 1.8471 56.4 0.1 
ThO, 1.609 1.34302 1.3325 2.0893 1.9834 1.9863 103.0 —2.9 

S 3.805 1.35917 1.33408 2.0985 1.9901 1.9934 105.1 —3.3 

S 0.395 1.33592 1.33329 2.1070 1.9956 2.0000 107.0 —4.4 

TiOz 0.373 1.33657 1.3325 2.8599 2.3518 2.4236 436.5 —71.8 

TiOz 7.577 1.41991 1.3325 3.0080 2.3878 2.4861 521.9 —98.3 








TABLE IV. Refractive indices at various sol concentrations. I. Si0:—sol. 








(nya) —np(5)) — ("p(9) —p(5)) 
103 103 






































¢ X10? ne Nm Np(8) Mp9) p(s) 
0.5815 1.33353 1.3325 1.5144 1.5097 1.5099 4.5 —0.2 
0.8626 1.33394 1.3325 1.5037 1.4995 1.4990 4.7 0.5 
1.0213 1.33421 1.3325 1.5042 1.5000 1.5000 4.2 0.0 
1.3155 1.33469 1.3325 1.5032 1.4991 1.4990 4.2 0.1 
1.6076 1.33514 1.3325 1.5008 1.4968 1.4966 4.2 0.2 
1.9460 1.33561 1.3325 1.4962 1.4924 1.4923 3.9 0.1 
2.5815 1.33634 1.3325 1.4846 1.4813 1.4813 3.3 0.0 
3.0751 1.33705 1.3325 1.4837 1.4805 1.4803 3.4 0.2 

TABLE V. Refractive indices at various sol concentrations. II. S—sol. 
(np) —M9)) = (Mpc) — (5) 
g S/100 cc ne nm p(s) Np(9) Np(5) x<103 x<103 
0.79 1.33592 1.33329 2.1070 1.9956 2.0000 107.0 —4.4 
1.59 1.33871 1.33338 2.1136 2.0001 2.0038 109.8 —3.7 
2.72 1.34256 1.33347 2.1108 1.9982 2.0022 108.6 —4.0 
3.81 1.34612 1.33364 2.0921 1.9856 1.9890 103.1 —3.4 
4.23 1.34755 1.33370 2.0918 1.9854 1.9887 103.1 —3.3 
4.76 1.34933 1.33375 2.0915 1.9852 1.9882 103.3 —3.0 
5.44 1.35171 1.33384 2.0949 1.9875 1.9908 104.1 —3.3 
6.34 1.35473 1.33394 2.0933 1.9865 1.9899 103.4 —3.4 
7.61 1.35917 1.33408 2.0985 1.9901 1.9934 105.1 —3.3 
TaBLE VI. Refractive indices in methyl-ethylketon-sols. 

pan ¢ X10? me Nm N p(s) N p(9) Mpls) anes ae ae aaa 
Acetylcellulose 1.047 1.37899 1.3778 1.4933 1.4914 1.4919 1.4 —0.5 
Acetylcellulose 7.277 1.38711 1.3778 1.5080 1.5056 1.5057 2.3 —0.1 
Nitrocellulose 2.746 1.37841 1.3741 1.5346 1.5309 1.5310 3.6 —0.1 








order to make the probable error negligible. This error will then often not be negligible. For ex- 


is particularly true in the case of colloids which ample, if water is the medium, the error will most 
are stable only in one medium, e.g., in an aqueous often be in the fourth or third decimal if the 
medium. A selection of the medium is excluded dispersed substance is a liquid or a solid re- 
also in those cases where the value of m, shall be spectively. 

used for obtaining information on the approxi- A survey on the variation of the error with 
mate degree of solvation (hydration) of the (#,—m,) is given in Table III for aqueous col- 
particles in a particular medium. The probable _loidal systems. Experimental data on my and m» 
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Fic. 1. (Wherever N appears in the figure m should be substituted.) 


as obtained by Freundlich and Rogowski" and by 
Odén” are used. It is assumed that the n,-values 
obtained by means of Eq. (5), m:s), are correct in 
a first approximation and that consequently the 
difference (ms) —mp s)) can be taken as an ap- 
proximate criterion on the error committed in 
applying Eq. (8) in the case of large (n,—n,») 
values. 

An inspection of Table II] shows that the log 
of the fractional discrepancy (8) —:5)) /M pcs) is 
a linear.function of the log of ( :s)—mm). This 
relationship is plotted in Fig. 1 which contains 
data of the Tables IV and V. From Fig. 1 the 


empirical function 


Np(8) 
Np(9) = (9) 
0.093 (p(s) — %m)!*° +1 ; 





is obtained. The use of this correction equation 
reduces the discrepancy considerably if (2:3) — %m) 
assumes large values. For small (15:3) —») values, 
on the other hand, the use of the correction 
equation is not advisable. 


F. Influence of the Colloid Concentration 


The influence of the colloid concentration is 
tested in Tables IV and V, the latter referring to 
a case where (3) —%m) is far too large to give 
reliable results by means of Eq. (8) alone. In 
either of the two systems, the calculated n, 
increases somewhat with decreasing colloid con- 


1H. Freundlich and F. Rogowski, Kolloid Beih. 37, 215 


(1933). 
2S. Odén, Zeits. f. physik. Chemie.80, 735 (1912). 


centration. The increase is the same whether 
Eq. (8) or the equation combination (9) or 
Eq. (5) are used. It is therefore advisable to 
calculate m, for infinite dilution. Differential 
measurements of (m;—m,,) as directly used in 
Eq. (8) are the only procedure which will allow 
one to obtain the most accurate m, for g—0. 


G. The Influence of the Absolute n,,-Values 


Careful refractometric measurements on non- 
aqueous colloidal solutions are scarce. A few data 
of Freundlich and Rogowski are used in Table VI. 
The comparison of these data with those of 
Table III seems to indicate that the absolute 
value of ,, has no influence upon the validity of 
Eq. (8). No influence was to be expected. 


Il. DETERMINATION OF REFRACTIVE INDICES 
FROM “ABSORPTION” MEASUREMENTS 


Equation (6b) can be transformed to 


Nn? +2f 
ny? = N»*———_, (10) 
Nm? — f 
where 


f Ao? T i ~?( T ) ~t/ T ) 
2 ——) ~ am\6x*»d =. 2s N6x*mc 


(m: the mass of a particle, p its density, c the 
concentration, in g cm~*). Since 


I/Ip=7-1, 

















the refractive index m, can be determined from 
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measurements of the apparent (conservative) 
absorption without any refractometric measure- 
ment being necessary, in the case of known mm, 
provided only that the particle weight (molecular 
weight), the density and the concentration of the 
particles are known. Since data on the two latter 
quantities must be known also if any of the 
preceding mixture rules are to be applied, the 
practical difference of Eq. (10) compared to the 
earlier ones is that the m; or (m;—m,) measure- 
ment is replaced by an absorption measurement. 
This practical difference defines the cases where 
the use of Eq. (10) is of particular practical 
interest (very turbid systems or very opaque 
particles). 

This alternate principle of m, determinations is 
restricted, of course, to application under the 
well-known conditions (particle size, particle 
shape, colloid concentration, and value of 
(m,—Mm)) for which the Rayleigh equation is 
valid. If all conditions are met except the condi- 
tion of spherical shape and of isotropy, Eq. (10) 
may still be applied, however, by introducing the 
Cabannes factor." This means that r must be 
replaced by 
1—(7/6)6 

1+0.55 


T 


where the depolarization factor 6=7,,/i, (7, is the 
intensity of the ‘Rayleigh component” which 
vibrates perpendicular to the plane of observa- 
tion and 7, is the intensity of the orthogonal 
component). 


APPENDIX—OPTICAL DETERMINATION 
OF DENSITIES 


The density of colloidal particles is generally 
determined by pyknometric methods or by 
methods of floating. In some instances, it may be 
of advantage to substitute optical methods for 
these methods. There are several possibilities 
which are implied in equations given in preceding 
sections of this paper. 

The transformation of Eq. (5) leads to 


Pp=Pm— 
Wm Ni—Nm 


Wp Np—Nm 
= pr— 








(11) 


13 J. Cabannes, Ann. de physique 15, 5 (1921); La 
. diffusion moléculaire de la lumiére (Paris, 1929). See also 
P. Debye, reference 7. 


where w is the weight in g. The use of this equa- 
tion differs from the use of the conventional 
Eq. (12) (case of exclusive pyknometric determi- 
nations of p») 


Wy 
(12) 





Pp 
W:/ pr - Wm/ Pm 


by the substitution of two gravimetric-volumetric 
data, e.g., of wm and p» by two refractometric 
data (m;—mm) and ny. Both Eqs. (11) and (12) 
should be valid on the assumption of additivity 


and of negligible volume contraction. 
More interesting appears the relation 


Nm( Ny? — Nm?) 


, (13) 
(mt —MNm) (Mp? +2Mm?) 





3 
Pp=-—C 
2 


which is obtained by transformation of Eq. (8). 
This equation makes it possible to determine 
densities exclusively from refractometric data, 
provided, of course, that the concentration is 
known. The application of Eq. (13) is subject to 
the restrictions imposed on the use of Eq. (8). 

If the particle weight (molecular weight) and 
the light scattering power of the dispersed sub- 
stance is known, then part of the refractometric 
data may be replaced by data on turbidity 
(absorption). This is possible by transforming 
Eq. (6b) to: 

243 Ny (Ny? — Nm?)? m 


—- (14) 


Pp ?» 


Not (my2+2ttm?)? 7 


2 Ny? —Nn? (6r'cm\! 
p= — at") . (14a) 
Ao? Mp? +20? T 





or 


Equations (13) and (14) have the disadvantage 
of requiring that m, is known. This will often be a 
deterrent to using them in the place of Eq. (12). 

This disadvantage can be avoided by using 
rather a transformation of Eq. (7) 


3ro4 T Cc 
oe (15) 


Pp= 
” 32s Nm? (Ny — Nm)? v 





which requires that the volume of the particles 
(molecules) is known. 





equa- 
‘ional 
ermi- 


(12) 


etric 
etric 

(12) 
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The boundary conditions of free vibration can be satisfied on the major surfaces of a plane- 
parallel plate if the displacement components are assumed to be products of trigonometric 
functions. In addition, the boundary conditions can be approximately satisfied on the minor 
surfaces when the plate is thin. The theory leads to a frequency equation 


v= 4(c/p)§[(m/2b)?+k(m/2a)*}? 


which has been found empirically to satisfy observations. The theoretical values of the constant 
k are 3.7 and 1.8 for the AT and BT quartz plates, respectively, while the observed values are 
3.9 and 1.7, respectively. 

















INTRODUCTION ; those quartz plates whose one major edge is 


REE elastic vibrations of anisotropic plates parallel to the crystallographic x-axis. Most of 

were treated extensively by W. Voigt! who the experimental work was done on this type of 
developed a general theory. Voigt’s theory is quartz plate. 
limited to the case where the wave-length is large Standing sine-waves along one major edge will 
as compared with the thickness and it must fail be first considered so that the boundary con- 
when the wave-length is of the order of mag- ditions on the major surface can be met rigor- 
nitude or smaller than the thickness. This latter ously. With these solutions, the boundary con- 
case was not interesting to Voigt; because of their ditions on the smaller lateral surfaces can be 
high frequency it was impractical to excite these Satisfied approximately. 
vibrations. By piezoelectric excitation, it has Standing Sine-Waves along the X-Direction 
been possible to excite vibrations of high fre- 
quency in thin plates. I. Koga? treated them 
theoretically by assuming the plate to be 
infinitely extended and gave a rigorous solution 
for this case. Experiments** showed, however, 
that while one resonant frequency and the cal- 
culated mode of motion were in good agreement, 
numerous other frequencies could be detected C15 = Cig = Cog = Cog = Cap = Cag = Cag = Cag = 0, 
which were not foreseen by the theory of the the cinemnateaia' ueletionn sna? 
infinite plate. The purpose of the present paper 
is to explain these deviations by presenting a 











We consider a rectangular plane-parallel plate 
of edges 2a, 2b, and 2c where 20 is the thickness. 
The origin of coordinates is in the center and 
the x-, y-, and z-directions are parallel to the 
edges 2a, 2b, and 2c, respectively. If, as we 
assume, 










—Xs =CiXetCiVytCiszetCisys 







theory of the thin plate whose lateral extension — Yy=CisXz2+Co2VytCostet+Couye 
is large but not infinite. —Z,=CisXet+CosVy Caste +Casye 
The consideration will be limited to the case (1) 





— Vye=CygXet Cosy tCoaBstCaay: 
—-X,= CosXetCoexy 





where the elastic constants C15, Cie, Cos, C26, Cas, 
C36, C45, and Cag are zero or negligibly small. Ac- 
cording to Voigt,’ the former is the case when —Xy= CoeXetCoery. 
the x-axis of the reference system (which we he boundary conditions for the free body are: 
assume parallel to one long edge) coincides with 

P g edge) X, cos (n, x) +X, cos (m, y) +X, cos (n, z) =0, 


a twofold axis of symmetry. This case includes ¥ y ¥ ao « 
—_—_—- cos (m, x cos (, . cos (m, z) =0, 
1W. Voigt, Lehrbuch der Kristallphysik (B. G. Teubner, P ( +% (m, 9) + (1, 8) (2) 










Leipzig, 1928). X,cos (n, x)+ Y, cos (n, y)+Z, cos (n, z) =0, 
21. Koga, Physics 3, 70 (1932). —_—_— ‘ 
$j. v. Aaa and P. J. Hart, Phys. Rev. 59, 85 (1941). *Voigt’s notations are used, so that compressional 
*'W. P. Mason, Bell Sys. Tech. J. 19, 74 (1940). stresses are positive. 
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which should be satisfied on all six faces. The obtain: 


dynamic equations are: ~X, steve +cupV¥ +00) 





























it cen Priva Poon Xexp [i(ax+py) ], 
dx Oy o8 — Vy =i(C:aU +c22p V+copW) 
Pa Phnsns Pond (3) xexp [i(ax+Py) ] 
dx dy az —Z,=1(C13aU +623) V-+Cs4pW) 
er a nn 
Ox Oy d8 — V,=i(¢yaU+cupV+cupW) 
where w=2zy is the angular frequency and 4, », Xexp [i(ax+py) ], 
and w are the components of the displacement 


vector. —X,=1(C5saW+ csp U+CseaV) 


We desire at first to satisfy Eqs. (3) in the Xexp [i(ax+py) ], 
whole body and Eggs. (2) on the major surfaces) — X,=i(cssaW+ceopU+cosaV) 
y=+b only. Furthermore, the displacement a 
components will be assumed to be independent Xexp [ilax+ py) } 
of z. Tentatively, the solution will be written Substituting into Eqs. (3), we get: 
li bination of plane w f the 7 ‘ ‘ . 
as a linear combination of plane waves of the Uh ~pelh hcal +eud + Veblen tod 


form 
+ Wap(cis+cs6) =0, 


u=U exp [i(ax+py) ], 
U . OY so all 2 2 
ee ee (4) ap(Ci2+Ce66) + V( — pw?+ Cosa? + Coop”) (6) 


W 2 op”) =0, 
w= W exp [i(ax+py)], _ Serna 
Uap (C1 +C56) + V(Csea? +Cosp”) 
where U, V, and W are constants. 


. ‘ W(— pw? +550? +C44p") = 0. 
If the strains are calculated in the usual way + W(— peat + coset’ + Cap?) 
from Eqs. (4) and substituted irito Eqs. (1), we The compatability of Eqs. (6) requires that 





CesP? +1107 — pw” ap(Ci2+Cee) ap(Cist+Cse) 
ap(Ci2+Cee) Coop* + Ceea* — pw” C5ea? +Cosp” =(). (7) 
ap(Cist+Css) C560? + Cop” Casp? +5500 — pw” 


Solving Eqs. (7) for p?, we obtain three roots p,’, p,., p3?, and from Eqs. (6) the corresponding ratios 
Ui: Vi:Wi; Us: V2: We and U3: V3:Ws, where the p,’s and the ratios U,: V;: W; are functions of 
pw” and a’. 
The form of this function can be seen if we substitute p/a=q and divide each row of the deter- 
minant (7) by a?: 
Ce6Q? +011 + pw?/ a? Q(Ci2+Coo) § BESS” 9(C1st+Cse) 
q(Ci2t+Cee) C229’ + C66 — pw*/a? Coot Cog? =(. (8) 
q(C14+Cs6) Cae t+ Coq’ CasQ? +C55— pw?/a? 





From this, we see that the roots p;7 will have the waves: 
form u=U; exp [i(ax+p,y) ] 
p2=a'f;(pw*/a’), (9) oa ee ae ee ee ee 
and ; 
and the ratios U;: V;:W; are only functions of u = U; exp [i(ax—piy) ]. 
pw*/a®. To each value p,? belong two complex eee er er 
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We form real linear combinations of these com- 


plex waves: 
u =U; cos ax sin piy, 


v= V; sin ax cos px, (10A) 
w= W; sin ax cos piy; 
” u= U; cos ax cos piy, 
v= V; sin ax sin py, (10B) 
w= W; sin ax sin py; 
and the final solution will be a linear combination 
of the three functions (10A) or (10B): 


u=cos ax:>, L;U; sin piy, 
v=sin ax->, L:V; cos piy, (11A) 
w=sin ax: >, L;W; cos py; 


nl u=cos ax: >, L;U; cos py, 
v=sin ax-)> L;V; sin pi, (11B) 
w=sin ax->, L;W; sin piy; 

where the U;--- and p; are known functions of 


pw* and a?, while the L; are to be determined by 
the boundary conditions. The strains are now: 


x,=—asin ax: >, L,U; sin px, 
Yy= —sin ax: >> p.LiV;i sin pry, 
Zy= —sin ax: >> p:LiW; sin piy, (12A) 
z,=a cos ax:>, L;W; cos piy, 
xy=cos ax: > L,(UipitaV,) cos px, 
for case A; or 
x,=—asin ax: > L;U; cos py, 
yy=sin ax:>> piLiV; cos pry, 
zy=sin ax->> piLiW; cos py, . (12B) 
z,=acos ax:>, L;W;sin px, 
xy=cos ax: >> Li(aVi—p,U;)) sin py, 
for case B. 


The boundary conditions require that the 


stresses 
Y,=Y,=X,=0 at y=+b. 


Substituting the strains, Eqs. (12), into Eqs. 
(1), we obtain 


Dd LileizaU ite piVitcap:W;) sin pib=0, 
Dd LilersaU i +cupiVitcupiW,) sin pib=0, (13A) 
Dd LilcespiU i+-CosaVitcCseaW;) cos pib=0, 


for case A; or 
dD Lil —Cr2aU + Co2piV; 
+ copiW;) cos pib=0, 
> Lil —cua +cupiVi 
+cupiW;) cos pib=0, 
DY Lil —cospiVitcesaV; 
+csseaW,) sin pb =0, 


(13B) 


for case B. 

The compatibility of Eqs. (13) demands that 
the corresponding determinant vanish. 

In view of Eq. (9), we can divide each Eq. 
(13) by @ and we then have only functions of 
pw*/a? in the parentheses. The transcendental 
functions sin p,b and cos p,b will have the form 


sin abf}(pw*/a’), cos abf(pw*/a’). 


Consequently, the determinant of Eqs. (13) will 
yield a transcendental equation which relates 
pw*/a® to ab. Solution of this equation yields the 
frequency as a function of the wave-length 1/a 
and of the product ab under the form 


pw” = a°g(ab). 


When this is done, Eqs. (8) yield the ratio 
L,:L2:L3 and the displacement components are 
known from Eqs. (11). 

While this formal solution. is impractical for 
numerical results, a numerical solution can be 
carried out by assuming some numerical value 
for pw*/a® and calculating the roots, q;, of Eq. 
(8) and the corresponding ratios U;: V;: W;. 

Then Eqs. (13) yield a transcendental equation 
for ab, the only unknown in (13). Thus the whole 
frequency spectrum can be determined. _ 

It is shown in Appendix I that Eqs. (7) and 
(13) include as a special case Timoshenko’s 
theory of flexural vibrations in an isotropic plate. 


Thickness Vibrations 


We shall discuss Eqs. (7) and (13) for the case 
where a is very small, i.e., the displacement is a 
very slowly variable function of x. Physically, 
this can be expected.to be true in a plate whose 
thickness is very small as compared to its lateral 
dimensions. If the displacement is independent 
of x as can be assumed in an infinitely extended 
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plate, a=0 and Eq. (7) becomes 
CesP? — pw” 0 0 

0 Coop? — pw? Cup? |=0. (14) 

0 Coup ® 


This determinant has the characteristic values 











Por? = pw*/Ces, Por? = pw*/C2’, pos*=pw*/cs’, (15) 
where 
Co,3' =3(CootCast[ (Co2—Cas)? +4024? }#). (16) 


I. Koga? has shown that the boundary conditions 
on the surface of an infinitely extended plane- 
parallel plate are satisfied when 


on panes ) = 
iti PT Oe 







This must also be a solution of Eqs. (13A) or 
(13B). Indeed, if 


Poi = (pw?/ces)t'=nx/2b (nm odd), (18) 


then we use the functions (11A). From the deter- 
minant (14), it is evident that the corresponding 
ratio is 







Ui: Vi: Wi =1:0:0. 
Because of Eq. (18) 
cos bob =0. 







Therefore, the Eqs. (13A) are satisfied by 
Li:L3:L3= 1:0:0. 








If, on the other hand, 
Poo=nx/2b or Pos=n2/2b (n odd), 


then we use the functions (10B). From Eqs. (6) 
and (14), it can be seen that in both cases 













Us = Us; =(). 
Again, 
cos Pood =0, 
or 
cos po3b =0; 
and (13B) is satisfied with 
| Ly =L;=0, 
or 
Ly =L, =0, 






respectively. Similar considerations show that for 
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n even, Koga’s solution satisfies Eqs. (13A or B) 
if a=0. 

We return to the case where a is not zero, but 
small. Equations (6) can be written as 


(A +a/pB+a°/p?C)U=pw*/p'U, (19) 


Cee 0 0 
O Coo Coa}, 
O Coa Cas 


0 CiotCes Cis t+Coe 
B= C12+Ces 0 0 , (20) 


where 


A 


Cua tCoe 0 0 
Ci 0 0 
C= 0 Ces C56 F- 

O Cse Cos 

Putting 
1/pB+a/p?C=D, (21) 

Eq. (19) becomes 

(A+aD)U = pw*/p’U. (22) 


Disregarding for the moment the fact that D 
still contains a and p, we obtain from perturba- 
tion theory the second-order solution of Eq. (22): 


DyiVor 








men " pw?/ Dor? — pu*/poi? wn 

pw pw* : Di? 

as are “ 7a*/ Po? — pu*/ Par (24) 

where 

Dir = (DU ox: Vor) = (Woe: {1/PB+a/p*?C} Uo) 
=1/p(Uox.- BU) +a/p?(Uox: CU) (25) 
=1/pByita/p?Cui. 


In Eq. (25), Uo, are the normalized eigenvectors 
of the secular equation (14) which satisfy the 
following relations: 

C22 Voot+CosWo2 = C2’ Vor, 

Cos Vor t+CusaWor= Ce’ Woe; 

C22 Vos+CosWos=Cs' Vos= —C3' Wor, 

C24 VostCasWos=Cs' Wos=Cs’ Vor. 


(26a) 


(26b) 


The third member of Eqs. (26b) follow from the 
orthogonality of the eigenvectors. Further, 


Vo = (1 ’ 0, 0), Vos = (0, Vo, Woz), 


7 
Uos=(0, — Woe, Vor), - 

















3) 


ut 


9) 


0) 














with 


Co2— C2’ \? + 
Voo=} 1+ = W;, 
C24 
(28) 


(C22—C2’) C22 — Ce’ 4 
Woe= — . 1+ = Vos; 
C24 C24 


and, of course, 


Voor? + Wo?’ ao Vos? + Wo;* =1. (29) 











The expansions (23) and (24) will fail if the 
matrix A is accidentally degenerate or quasi- 
degenerate, i.e., if two of the values 

Ces, C2’, Cs 
are equal or nearly equal. In this case, the de- 
generacy should be first removed by the con- 
ventional methods. 

It can be seen that B,, =(Uo,.-BUo,) vanishes. 
Indeed, the vectors Uo, have either no y- and 
z-components or no x-component. Application of 
the matrix B to these vectors produces a new 
vector which has no x-component or no y- and 
z-components, respectively, so that the scalar 
product of vector and transformed vector 
vanishes in either case. Therefore, by Eq. (25): 


Dix=a/p?(CUon-Uor) =a/p?Crr. (30) 
According to Eq. (25): 
Di? =1/p?(Wox- BU)? +--+ =1/p?B? 
(k#1), (31) 


where the dots stand for terms of higher order. 
Equations (23) and (24) become 
, ByVo 
U.=Uorta/por DU - 7 » (32) 
pw*/ Por — pw*/ por’ 


pw/ pi? = pw*/ por? +a?/ por Ke (33) 
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2 (Cis t+C56)C2a(Ci2+Ces) — (C1a+C56)?(Co2 — Cos) — (Ci2-+Co6)*(Cas — Cos) 





PLATES 


where 


Bi? 





kKe=Crr+ >’ —, 
pw/ Pox? — pw*/ por 

In Eqs. (32) and (33) we have, on the right- 
hand side, replaced p,2 by po,” because a cor- 
rection would only add terms of higher order. 
As shown above, if a=0, 


[Ces 

pw? /B? = + ce’ 
/ 

C3. 


If a@ is small, the expansion 


[ cos 
pu? /B8? = (1+e)+ ce’ (34) 


lcs 


can be assumed. ¢€ must be determined as a 
function of a. Combining Eqs. (33), (34), and 
(15), we get 
be pw B?(1+€)c 
; CHarKe/pow Ct+a°cK:/ pw* 
 BL+0c 
~ C+0°,/B*(1+6) 


where c stands for either one of the three values 
given in Eq. (34), and the notation 


a/B=r, (36) 








8(1+e—r°x,/c), (35) 


is used. 

In the following, we shall treat explicitly the 
shear vibration whose zero-order approximation is 
given by Eq. (18). If 


pw*/B? = (1+ €)cee, (37) 
then by Eq. (35) 
pr=B(1+e—7'x1/Ces). (38) 


xk, is given by Eq. (33). In appendix II it is shown 
that explicitly 






(39) 





Ki =Ci1 


(C22 — Ces) (Cas — Ces) —Co? 


Substituting the expansions (32) and (33) with (27) and (17) into (13A), we have 





B 
Licee(Poit ° 4.) cos pib+Lra( co r ~ 
Co — C66 
B 
+Lsa (cu ; > 
C3 — Cee 





+ Cee VostCseWost : - ‘) cos (posb+:-- 





+66Vo2+cseWoet : - -) cos (po2b+---) 
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Vi ByW. 
Lra( cuter > =e +eu > = “) sin Poib+L2(Co4Vo2+CssWoet - - 
2,3 Cos — Ci’ 2,3 Ce6— c,’ 
X (poo t:-+) sin (poob+--+)+L3(CosVost+CsaWos) (post - - 
Sule | By.W. 
Lya( cater % Cu > “) sin £1b-+Le2(C22 Vo2+CosWoet : - -) 
2,3 “ag 2,3 Cee— Cr 
X(poet+-+-) sin (poob+ ---)+L3(Co2Vos+CosaWost+:--) (post: : 


The dots stand for terms of higher order which can be omitted. 
The terms of Eqs. (40) will be simplified one by one: By Eq. (38) 


The upper sign refers to n=1, 5, 


cos p1b=cos =( +4[e—r'x1/Ces |) EF 


pi:28(1+43[Le—r?x1/Co6]), 


-) sin (posb+ coe 


-) sin (posb+ oe 





a (40) 













9-.-- and the lower to »=3, 7, 11-- 


appear with the first terms of the second and third Eqs. (40) are 


The two signs refer again to the two cases mentioned above. Smaller terms are omitted here as they 


are in addition to unity. By Eqs. ( 


By.= (Uni: 
Bi3= (Uo1- BU os) = — (C12 +66) Wo2t+ (C1a+C56) Vor; 


Bis 
on. ow Voo+CseWo2= Von cus) + Wea ee 


Co’ —Coe 





Cee 


C3 — Cee Cc: = 








Bis 
+e) roa + Css Wos = Cos——_—— —CosWo2tCseVo2 


sin p,bSsin nx/2= +1. 


20) and (27) 
BUoz2) = (€12-+€66) Voot+ (C14 +56) Woe, 


Co — Cee C2 — Cee 





C2 — Cee 


Cee 





mm 
7 1/08). 


-. The factors sin ~,b which 


) 


1 
=— [ Voocee(C12+C2’) + Wo2(CrsCoe+C2’Cse) |; 


(42) 









(43) 







(44) 






(45) 






































= —[ Voo(CesCia+Cs’Cs6) — Woetes(Ci2+cs’) ]. (46) 
C3’ — Cos 
In view of Eqs. (26a) and (26b): 
Bi2Voe 4, Pu Vos BuWes | BisWos Co'BypWoe C3’ BisWos 
Cia tC = r + - Ci4. (47) 
ocad”, Cos— Cs a —C; Cee — C2 Cee — Cz 
The second and third terms of the second Eq. (40) become in view of (26a) and (26b) 
pone L262’ Po2Wo2 sin poob+Laes' pos Voz sin posd. (48) 
Bie = ia.. Vos wi BisWos Co'By2Vo2  C3'BisVos 
Cia +C22 >; Hee t - . (49) 
C66 — 66— C3" C66 — mie C66 — C2 Ces — C3 


The second and third terms of the third Eq. (40) become in view of (26a) and (26b): 


L2l2 
According to Eqs. (15) and (18): 


pu=B, 


‘Por Sin Porb — Les’ Woofos sin posd. 


Porz=B(Cee/C2’)*, pos=B(Cee/Cs’)?. 
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Equations (40) read now: 


a 


Co’ — Cos 


nn n 
FLiceb——(¢ —rx/Ces) +L [ Vo2Cee(Ci2+C2’) + Wo2(CraCoo+Co' Css) | cos 5 (cos/ex') 


a n 
+L; : [ Vo2(CisCoe+Cs Coe) — Wo2cee(Cr2 +c;') ] COs 5 (ew/ea) =0, 


C3 — Cee 


Co’ BizWo2 C3’ BisWos 
+ 


, 





nT 
+1ra(cut ) L2€2'B(C66/C2’)' Woe sin p (ces/ex')! 


Cos — C2" Cea — C3 


(50) 
ne 
+ L3C3'B(Ce6/cs')* Vo2 sin py (eea/en)? =0, 


C2'Bi2Vo2 Cs’ Bis Vos 





Cos — Co’ Cos—Ca 


n 
+lia (c+ ) + L2C2'B(Ce6/C2’)* Vo2 sin 5 (ces/e!) 


nn 
— LCs’ Wo2(Coe/cs’)' sin > (ces/ea')! =0. 


Consequently, the determinant 


nn a nn , 
¥ Cos8—(€—1?x/Cos) - [ Vo2ces(Ci2+C2’) a cos —(C¢66/C2 ) 
4 C2 — Coe 





- [ Vo2(C1sCoe+Cs’ Coe) 
+ Wo2(CrsCoe+02'c )] cos —(c /c2')3 acl 
02(C1sCo6 +Co' Coe 7 (ae/es — Worcedlernter')] 


C2’ By2Woe 4 , na . _ nr 
ta(cut C2'B(Ce6/C2’)*Woe sin > (esalen')* C3'B(Ce6/Cs’)*Vo2 sin > (oases)? 
66— C2 


oa) 
+— — 
C66 — C3 


C2’ By2Vo2 ‘ NIV. of nn 1. #N , 2 ne oe ay 
tal ¢12.-+—— _ C2'B(Ce6/C2')* Voz sin > ew C2’)? —C3'BWo2(Cee/Cs ) sin y (eases ) 
C66 — C2 


C3'BisVos 
— 
Cea — C3 
must vanish. The determinant (51) is expanded in terms of the elements of the first row. The first 


co-factor is 


ot _ nr Ore, ’ 7 
Ai =Co6(Co'C3 )is? sin y (eaeles )! sin > (easles )( — Wo2? — Vo2") 








nT nT 
= —C66(Co’c;')§8* sin > (ees/ea')! sin > (eases (52) 


in view of Eq. (29). The second co-factor is: 


‘ ‘ Co BizaWoe C3’ BisWos 
Ai2 = tacs'B(Co6/cs’)| { Crat — - 02 


C66 — Co Cee— C3 


'BiasV. "Bis V. 
+(o+" 12 a6 13 ) Vea] sin (cu/e') 


Cos— Co" Cos— Cs 








cB c;'B 
= sacya(eu/cr')} cases (Wo2?+ Vo2*) + —- (Wo2Wost+ Voe Vos) 


C68 — C2" Cos— Cs" 


nr 
+¢14.Wo2tC12 Ve sin p (eae/es). (53) 
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The second term in the bracket vanishes because the quantity in parenthesis represents the scalar 
product of two orthogonal eigenvectors. We substitute the value of By. from Eqs. (44) into (53): 


Co (Cia + Co (C te n 
dum tap (cwcr)'| Vou( cut BE) 4 +Wa(et = : A = See |i 5 (coe/es') 
C2 


Cee 66 — C2" 











(Coscs’ ) 


mnt Y ‘o2Cee(Ci2+Ce’ )A Wo2(CrsCoo +02’ C56) | sin ~ (coe/Cs’ yh (54) 


C6 6— C2’ 








= +a 






The third co-factor is: 





Co'Bi2Woe 3 Bis 2) 
+ 


Cos — Co" Coa —C3" 


Co’ BizVoo C3’ BisW. nr 
“ee 2PizVor C3 Dis *) | si 5 (ces ca!) 


C68 — Co’ C66 — C3 


Ai3s= sacs (cu/ed) Voe (cut 























C3 Bis 
= saBer cu/cr)| -( Vo2? + Woe") iw ~(WoeV 02 — Wo2Vo2) 
Ces —C3' C6 6— Co’ 
‘ _ nT bi 
+ Voetius— Woecr| sin > ee C2) 
(C66/C2')! , , _ ’ ce 
= +afco nf Vo2(Cs’Cse+CraCes) — Wo2ces(Ci2+cs’) |] sin —(Co6/Cc2’)', (55) 
Cee C3 2 
because Vo2= Wo3 and Wo2= — Vo3. The determinant Eqs. (51) can now be written in view of Eqs. 





(52), (54), and (55): 







Tn ‘ i , un - , nT ‘ 
se Cee'B!-(e—1?a1/Cos) (C2 C3 )} sin > (eaes Ce )} sin > ee C3 ) 






3 nr nt 
(ceacs’)*L V o2€66(C12 +2’) + Wo2(CrsCoe-+C2'Cse) |? sin 9 ees c;')} cos > c;')} 





(c2’ — C66 





9 


a nr . nt 
Face ert Vo2(C1sCes-+C3' C56) — WoeCes(Ci2t+cs’) |? sin 9 ee C2’) cos Fhe /c2')'=0. 
C3 — C66 








Solving for e, we obtain e=r°k,, where 






nr 
cot —(C¢66/C2’)! 
1 4 


kn=—| mt ; —[ Voste(Ci2 +2’) + Woa(crscos+ce' Coe) }? 
n(Ces)! (Ce (Ces — 3 
















nT ‘ 
cot —(C6e/C3 )! 





[ Voo(CrsCost+Cs' Css) — Worcee(Ci2t+cs’) . (56) 





(cs’)*(Ces—Cs')* 





If ci4 and Cse are very small: 





4 wn 
, cot —(Ce6/Cc2’)! 
b 1 (C12+Ce6)?(C4s — Cos) Coe?(Ciz tC )? 4 2 

g=@—j Cu 


C66 (C22 — Coe) (Cas — Coe) — Coa” i+ (=) nm (Cee)? (c2’)4(Ce6— C2")? 


C24 





















n f / 
cot ? (Ce6/Cs')3 


+(4—*) Cos'(C12+es')? eee a | 
C22 — C2 ) nm(Cceos)! (cs! yi Ces—C3 )” 








C24 











C24 






W 
Vi 





lar 


54) 


55) 


qs. 


56) 


57) 
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The fourth term vanishes with ¢24. Indeed, by Eq. (26a), Wo2=0 when ¢c2,=0. Physically, the con- 
stant Co, guarantees that there shall be no purely longitudinal plane waves propagated along the 
y-direction, and that the displacement vectors Uo2 and Ups form angles with the wave vector. The 
fourth term in Eq. (57) is due to this fact. If co, vanishes, Eq. (57) becomes 
















nn 
Aces! (Ci2+C22)? cot —(Ce6/C22)! 
1 (c12+Co6)? 66 ( 12 22) 2 ( 66/ 22) 


k,=— Cu-— . + ’ (58) 
C66 Co2 — Cee NWC224(Coe on C22)? 

















because then cs’ =C22 and Vo2=1. In the isotropic case, 









C11 =Cop=2ut+A, Cis=A, Coo=H, 






while all other constants vanish, and Eq. (58) becomes 


j 
kp, -1+—(— ~~). (59) 
2u+ar 2u+aA 


We now return to the general case: The ratio L;:L:2:L; is equal to the ratio of the corresponding 
co-factors: 


L;: Le: L3= Ai: A12: Ais 
















vail om r V 02€66(Ci2+Ce’) + Wo2(Cr1sCoe-+Ce'Cse) | 
(C2 — Ces) (C2'Coe)! | 











_ nr : 
sin > (css/ee ) 





[ r Vo2(Cs'Cse-+C14C66) — Wo2ee6(Ci2t+Cs ) 
8 deen — . (60) 
(cs — Ces) (Cescs’)* 








re. 
| sin ? C66/C3 





If we substitute these values into Eq. (11A) and take into account Eqs. (32), (33), and (35), we 
obtain equations of the form 





u=cos ax sin By[1+}r?(k—Ki/Cee) |, 
v=rsin ax[A, cos By+A2 cos By(Ce6/C2’)'+As3 cos By(Ces/cs')* ], (61) 
w=rsin ax[A; cos By +A, cos By(Ce6/C2’)'+As5 cos By(Cos/Cs’)*], 







where the constants A;---As5 are only functions of the elastic constants. We do not need the explicit 
values, but it is simple to write them. 
In the isotropic case, L3, Wm, and Woz vanish. Therefore, we have 


r u(A+u) 
8 (_» ) (A+u)[u(A+ 2x) } 
A+2p 










w=0, Au: 412:413=1: + 











; 
cot =( | (62) 
2u+aA 2 \2u+a 


| w(*-) 
cos 
P \ 4 2u 


4 
v=r sin on —cos Ay-+2( ) ’ 
A+2y _ nT m ) 
sin — 
' 2\A+2u/ | 


8 
u=COSs ax sin ay +r—( 
nT 
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The difference between Eqs. (61) and the solution for the infinite plate is: 
(1) Instead of u=sin By we have a slow variation of u with x. 







(2) The derivative du/dy is not exactly zero at y= +b, because fp, differs by a small quantity 


(proportional to r*) from B=nx/2b. 


(3) While for the infinite plate » and w vanish, we find small quantities proportional to r in the 
case of the finite plate. v and w do not in general vanish anywhere. 


Approximate Solution for the Free Plate 


The Eqs. (1) and (13A or B) are rigorously 
valid for sinusoidal standing waves in a plane- 
parallel plate. Consequently, the frequencies (56) 
and the modes (61) are correct solutions for this 
case, if higher orders of r are disregarded. For a 
free limited plate, however, the boundary condi- 
tions in Eqs. (2) should be satisfied on the lateral 
surfaces. The boundary conditions in Eqs. (2) 
on the faces x = +a require that 


dv ow 
—X; meus ben teu 
dy 
=asin aa: F;(y)=0, 
—X, =e tou(— + 
oy ax 


(63) 
= £6 cos aa: F2(y) =0, 


sa ov 
—X,=Css— pete 


joo cos aa: F;(y) =0, 


where F;, F2, and F; are functions of y which, 
in general, are of order of unity. It is evident 
from Eqs. (63) that it is impossible to satisfy 
rigorously the boundary conditions by a choice 
of a. If the second and third Eqs. (63) are to be 
satisfied, cos aa must vanish, but then sin aa 
cannot vanish simultaneously, and therefore the 
first Eq. (63) will not be satisfied. However, if 
we put 

cosaa=0; a=mr/2a (m=1,3,5---) (64) 


the traction X, is proportional to 1/a and tends 
toward zero as the lateral dimensions are in- 
creased, while the other traction components 
vanish completely. We shall consider Eqs. (64) 
as the definition of a. By moving the origin of 
the reference system along the x-axis, we can 
write instead of Eqs. (61): 

u=sin ax sin yB-[1+47r?(k—K/ces) |, 

v=rcosax:[--- ], (61a) 

w=rcosax:[--- ]. 














The boundary conditions Eqs. (2) are now: 


—X,=a cos aa: F,(y) =0, 
—X,=—8 sin aa- F,(y) =0, (63a) 
—X,=-—8 sin aa: F;(y) = 


To satisfy the second and third Eqs. (63a), we 
have to assume 


sinaa=0, a=mr/2a (m=2,4,6-:-). (64a) 


Taking the even and odd modes of (64) and 
(64a) together, we note that m can now take on 
any integer value. 

If the extension in the z-direction, i.e., 2c, is 
not infinite, the boundary conditions 


Ou ov ow 
—Z2=Cis—+623—+C34- =0, 
Ox oy oy 
Ou ov ow 
one - — =0, (65) 
dy 


ow Ou dv 
—X, =Cos tees —-+— } =0. 
Oy Ox 


on the faces z= +c cannot be satisfied. However, 
if 
C13 = C14 = Cog = Cog = C34 = Coe =O, 

(e.g., in the isotropic case) w vanishes and the 
Eqs. (65) are satisfied exactly. Fortunately, these 
constants are often small so that we can expect 
a fair agreement between theory and experiment 
despite Eq. (65). 

We shall consider Eq. (61) and (61a) with the 
definitions (64) and (64a) of aas the approximate 
solution of the problem for the free thin plate. 


Comparison with Experiment 


Combining Eqs. (37), (56), (64), and (64a), we 


obtain: 
pw = BC e6(1+77k,) 


n\2 m \27 
r=eu/o) (=) +ba(™) | (66) 


Observations on thin quartz plates of the so- 
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called Ye type were carried out by Sykes.® In 
this case, one major edge is parallel to the 
twofold x-axis, so that comparison with the 
present theory is possible. Sykes found an 
empirical formula for the shear type vibrations: 


r= Heu/o) () *) (FZ) Jen 


When p=1, this empirical formula agrees with 
the theoretical ‘Eq. (66). The numerical value of 
the empirical constant & is not given explicitly, 
but it can be calculated from Sykes’ diagrams. 
Expansion of Eq. (66) yields, when »=1 and 
a>b: 





1 k PP 
Vn = A(coe/o)'-(1 a ’ 
and 
— k(b/a)?(m?—1 
oni faye) e44(/a) (mn? — 1); 
V1 2[1+4k(b/a)? | 


2 (: 2 tm — V1 
— =) : 
n?*—1\b Vi 


From Sykes’ Fig. 6.15, we find, for a/b=32: 
vi =1662, v3=1688, »,;=1736 


or 








which gives k=4.0 and 3.8, respectively. The 
average value is k=3.9 for the A7-cut plate. 
From Sykes’ Fig. 6.17, we find, for a/b=48: 


vy, = 2552, v; = 2600, vg = 2628, v1, = 2658 


which gives k=1.8, 1.7, and 1.6, respectively. 
The average value is k= 1.7 for the B7-cut plate. 
With the elastic constants of quartz measured 
by Mason® and the transformation formulas 
listed by Sykes,’ we obtain the following con- 
stants for the plates in units of 10" dyn/cm? 
(Table I). From Eqs. (56), (39), (16), and (28), 
we calculate the following values for k(n =1): 


AT BT 
Robs 3.9 1.7 
Reale 3.7 1.8 


The agreement is satisfactory. Theoretically, k 
should be independent of m, but depend on n. 
Sykes does not state whether this is the em- 
pirical meaning of k. However, comparison of 


*R. A. Sykes, Bell. Sys. Tech. J. 23, 52 (1944). 
*W. P. Mason, Bell Sys. Tech. J. 22, 178 (1943). 


VIBRATIONS OF THIN PLATES 

















TABLE I. 
AT BT AT BT 
Cu 86.05 86.05 C12 — 10.35 26.20 
C22 129.86 97.77 Cis 25.85 — 10.70 
Cus 39.06 40.85 Cu 3.55 0.13 
Cee 29.34 68.91 C4 —5.82 13.11 


Cue — 2.46 6.45 








Sykes’ Figs. 6.15 and 6.16 referring to the cases 
n=1 and n=3 seems to show that k is definitely 
lower in the latter case. 

It appears from the empirical Eq. (67) that 
for the case p=1, the frequency is not affected 
by the dimension c. This confirms our assumption 
that for the type of vibration under considera- 
tion, the boundary conditions (65) can be disre- 
garded. 


In the lowest mode (m=1), where 
a=n/2a, 


the displacement u is maximum in the center 
and decreases monotonically toward the edge 
x=-+ta. The displacement components v and w 
are small, but they have their maximum value 
at the edge. In a lycopodium powder experiment, 
we should expect the powder to be most vigor- 
ously moved at the center and only slightly 
at the edge, while no exact nodes appear. This 
is in agreement with experimental evidence.’ 
The clamping of a crystal on the corners does not 
modify considerably its mode, as can be expected. 
In the mode corresponding to m= 2, i.e., 


a=n/a, 


the mode is given by Eqs. (61a). The displace- 
ment component u has the form 


u=f(y) sin rx/a, 


which vanishes in the center while it has opposite 
signs in the left- and right-hand half of the plate. 
Observations by Sykes® (third figure in Fig. 6.5) 
are in agreement with this result. 

Experiments show that a freely vibrating 
crystal of the kind considered here, if mounted 
between metal plates parallel to its major faces, 
excites ultrasonic vibrations of the air between 
the surface and the plate. This cannot be ex- 
plained by the theory of the infinite plate 


7G. M. Thurston, U. S, Patent 1,883,111, October 18, 
1932. 
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because a shearing motion alone cannot cause 
compression or dilatation. The present theory 
shows that the component 7 is not zero and thus 
explains the existence of air waves. 

Sykes’ observations show that Eq. (67) ceases 
to be valid at certain ratios a:b and c:b. With 
these ratios, two nearby frequencies appear, one 
lower and one higher than that given by Eq. (67). 
Sykes explains this fact qualitatively by coupling 
of the thickness vibration with other types of 
modes. We shall discuss the phenomenon of 
coupling from the theoretical viewpoint. 

It was shown in a previous paper® how it is 
possible to represent the solution of an elastic 
vibration problem as a linear combination of 
“zero-order” or uncoupled modes. We ean con- 
sider the mode given by Eq. (61) as a zero-order 
mode. It differs from an exact solution in that 
it does not exactly comply with the boundary 
conditions on the lateral faces. Other zero-order 
modes can be obtained by different solutions of 
Eqs. (7) and (13). We have only considered 
those solutions where a is small. It can be shown, 
however, that other types of solutions corre- 
sponding to about the same frequency 


= i. 
v= (Cee/p) 4b 
can be obtained when a is of the order of mag- 
nitude of p. One type corresponds to Timo- 
shenko’s flexural waves treated in Appendix I. 
These other modes will be equally zero-order 
modes because they will satisfy the boundary 
conditions on the lateral faces only approxi- 
mately. Their frequency depends mainly on the 
dimensions a and c.° Consequently, a coupling 
in the sense defined in a previous paper® will 
exist between the thickness mode Eq. (61) and 
the “‘lateral’’ modes. As the coupling is caused 
by the stresses X, (on x= +a) and Z,, Y., and 
X, (on z= +0), which are all small, the coupling 
itself will be small. As shown by the form of the 
secular determinant in the same paper, a weak 
coupling has a noticeable effect only when the 
frequencies of two zero-order modes are nearly 
equal. This explains why the exceptional regions, 
where the uncoupled frequency becomes incor- 
. rect, occur at the intersection of the curve 


8H. Ekstein, Phys. Rev. 66, 108 (1944). 
* The proof is omitted here. 
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representing the uncoupled thickness shear mode 
and one representing an uncoupled “lateral’’ 
mode. 


APPENDIX I 


It will be shown that Eqs. (7) and (13) include, 
as a special case, Timoshenko’s theory of flexural 
waves in an isotropic plate.” In the isotropic case, 


C1 = Co =2u+A, Cu=Cos=Coo=mM, Cr=d, (1.1) 


and all other constants appearing in Eqs. (7) 
and (13) are zero. The determinant (7) reduces to 


up? +(A+2u)a®— pw? = ap(A+u) 
ap(A\+u) (A+ 2) p?+ na? — pw” 
with W,= W,=0. The third solution: 


up? + po? — pw* =0, 


=0, (1.2) 


and 

U3= V3;=0, 
can be disregarded for this discussion, as will be 
seen. The two solutions of (1.2) are: 


2 


nn} 1.3 
pi \42p a (1.3) 


p? os pw? /p —<@’. 





and 


If we substitute this into the first two Eqs. (6), 
the ratios 


U,;: Vi=a: pi, U2: V2=—peia (1.4) 


are found. We choose the case (A) corresponding 
to Eqs. (10A) and (13A). Equations (13A) reduce 
to 


> LildAaU;+(A+2u)p.V;] sin pb =0, 
> LippiW; sin pib =(), (1.5) 
¥ Lin(piUitaV,) cos pid =0. 


The second Eq. (1.5) can be satisfied identically 
by 
L;=0, 


because W,;=W:=0. We consider frequencies 
low enough so that: 


pw /p <a? 
and all the more: 
pw?/(A+2y) <a’. 
1S. P. Timoshenko, Phil. Mag. [6] 43, 125 (1922). 
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According to (1.3), 1: and p2 then are both Egs. (1.5) become, in view of (1.3) and (1.4): 











imaginary. If we put Li[\a? —m?(A+2y)] sinh mb] 
: pw” 4 +iL22una sinh mb=0, 
m=al 1 —-—————__ ], ; (1.7) 
| = a®(A+2yp) iL,2ma cosh mb+ L2(n?+ a?) cosh nb=0. 





» . (1.6) By eliminating ZL; and Lz from (1.7), we obtain 
5 

pw . o 66 . ” 

n  nencaes : Timoshenko’s ‘‘frequency equation’”’: 

aru 





I 
a 





4ua’mn tanh nb 
= (a?+-mn*)[ (A+2u)m?—da*®] tanh mb. (1.8) 






so that m and » are positive real numbers, the 












APPENDIX II 





Instead of using Eq. (33) for the explicit value of «1, it is easier to find it by direct expansion of 
the determinant (8). If 






p/a=q. (1/q<1) 


we divide the second and third rows and columns by q: 







g’( — pw*/ p? + Cee+C1:/¢") Ciot Cos Cia tCoe 
Cio t+Ceoe — pw*/p?+C22+Ce6/¢" Cost C5e/Q" =0. (2.1) 
Cis t+Coe CortCse/Q° — pw? / Pp? + Cys tC55/¢ 





In view of Eqs. (33) and (15): 





pw*/ pi? = Ceg+ aki / Por’, 


so that we can neglect the term with a*® everywhere except in the first element of the first column of 
(2.1), finally omitting the terms c;,/g’ where they are in addition to cy, we obtain: 


q’(— pw*/pi?+costei/g*) Crates CistCre 


Cio t+Coe Co2 — Cee Co «=| = 0. (2.2) 









Cit Coe C24 Css — Cos 





We expand (2.2) and get: 


q*( a pw” /p1?+-Cos+C11/Q") [ (Cor — Cos) (C44 — Coe) — Cx? | ae (C12 + C66) [ (C12-+Cos) (Cas — Ces) - Coa(Ciat+Coe) | 
+ (Cis+€56)[ (Ci2+066)Cos— (Cor —Cos) (Crates) ]=0. (2.3) 






Solving (2.3) for pw*/p,’: 





ps!/t=cuot1/¢ cn 





4 (Cia +56) {C24(Ci2+€66) — (C22 — Coe) (C14 +56) | — (C12 +Ce6) { (Ci2+€66) (Cas — C06) — Coa(Cra+Coe) } 


(Coo — Cos) (Cas — Coe) — Coa” 





| (2.4) 





In the expression: 1/g?=a*/p,? we can, in view of Eq. (15), substitute o,* for p:;? without adding 
correction terms of the second order. Comparing (2.4) with (33), we find that: 


2C24(Cis Cos) (C12+Cos) = (C14 +C56)?(Co2— Coe) — (Ci2+€o66)?(Cas — Cos) 


(C22 — Cee) (C44 — Coe) —C27 





(39) 





Ki=Cu 
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The effect of the hydrogen content of electrolytic manganese on the average lattice parameter 
was determined. The parameter varies approximately linearly with the hydrogen content, the 
enlargement in parameter being 0.0003 percent per cc of gas per 100 grams of metal. The 
crystals, composing electrolytic manganese metal, are generally not uniform in parameter but 
show deviations up to +0.2 or +0.3 percent. The deviations in parameter between the unit cells 
of alpha-manganese in a sample of metal are not proportional to the amount of gas in the metal. 
The effect of various degassing treatments on the parameter and diffraction patterns was 
determined. The parameters are not effected unless the metal is heated over 300°C. At 500°C 
and over, the parameter returns to its normal value and the diffraction lines are sharp. The 
broadening of the diffraction lines in the patterns for the metal containing gas are attributed to 
non-uniform lattices in the metal due to non-uniform gas distribution, and the results indicate 
that manganese may be saturated with hydrogen at about 695 cc per 100 grams of metal or 0.06 
percent. Manganese metal with normal lattice parameters and uniform structure may be ob- 
tained from electrolytic metal by removing the gas and heating it to 500°C for 1 hour. If the 
metal is heated to 500°C or more during the degassing treatment, subsequent heating is not 
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necessary. 








INTRODUCTION 


URING the rather extensive studies by the 
Bureau of Mines on electrolytic manganese 

and its alloys, a question has frequently arisen 
regarding the correct lattice parameter for alpha- 
manganese. The average parameter as deter- 
mined from several samples of electrolytic man- 
ganese metal was 8.904A. Moreover, the alpha- 
manganese phase in the Mn-Cu alloys‘ gave a 
parameter of 8.906A, so it was concluded that the 
alpha-manganese phase in the alloys was essen- 
tially pure manganese. The published® parameter 
for alpha-manganese was only 8.894A, but since 
the electrolytic metal was more than 99.9 percent 
pure and the impurities were distributed among 
a large number of elements, it was unlikely that 
the discrepancy could be due to impurities in the 
metal. There was no satisfactory explanation for 
the discrepancy until the hydrogen content of the 
electrolytic manganese was investigated. In this 


1 Published by permission of the Director, Bureau of 
Mines, U. S. Department of the Interior. 

2 Physicist, Western Region, Bureau of Mines. 

3 Metallurgical Engineer, Western Region, Bureau of 
Mines. 

*R. S. Dean, J. R. Long, T. R. Graham, E. V. Potter, 
and E. T. Hayes, ‘“‘The Copper-Manganese Equilibrium 


System,” Trans. A.S.M.E. 34, 443 (1945). 

5 Ralph W. G. Wyckoff, The Structure of Crystals, second 
edition and second edition supplement; Handbook of 
Chemistry and Physics, twenty-seventh edition, 1943-1944. 
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study some vacuum-distilled metal was made, 
and its parameter was found to be 8.8948A. This 
result checked the published value quite well, and 
the lines in the pattern were much sharper than 
those obtained from electrolytic manganese 
metal. Obviously, some impurity in the electro- 
lytic metal was causing the enlargement in 
parameter, but the amount of each of the known 
impurities was small when considered on a per- 
centage basis. The hydrogen content was quite 
high on a volume basis, however, being 250 cc per 
100 grams of metal,® so an investigation was 
made to determine what effect hydrogen might 
have on the parameters. 


TABLE I. Variation of average parameter with gas content. 














Hydrogen 
content in 
cc per 100 
Parameter grams 
Material in A of Mn 
Vacuum-distilled manganese 8.8948 +0.001 0 
Electrolytic manganese trans- 8.898 +0.005 95+13 
formed from gamma-manga- 
nese to alpha-manganese ; 
Regular electrolytic manganese 8.904 +0.004 250+25 
Electrolytic manganese, hydro- 8.910 +0.005 615+50 


gen content high 








‘ E. V. Potter, E. T. Hayes, and H. C. Lukens, “‘Hydro- 
= Content of Electrolytic Manganese and Its Removal,” 

m. Inst. of Mining and Metallurgical Eng.—Metals Tech- 
nology (June 1945), Technical Publication 1809. 
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CHANGE IN AVERAGE LATTICE PARAMETER WITH 
HYDROGEN CONTENT 


The parameters of a number of samples of 
regular electrolytic manganese were determined 
by the Debye-Scherrer method and the hydrogen 
content of the material was determined by the 
method outlined in an earlier publication.* Most 
of these samples had about the same hydrogen 
content ; one sample was found with an unusually 
high hydrogen content, and some samples of 
manganese, deposited as flexible gamma-manga- 
nese, had unusually low hydrogen content. The 
gamma-manganese changed gradually at room 
temperature to alpha-manganese, and the pa- 
rameters were obtained for the alpha-form. The 
results of these parameter and hydrogen-content 
determinations are shown in Table I along with 
corresponding values for vacuum-distilled man- 
ganese. 

These data are plotted in Fig. 1 and show a 
linear relation between parameter and hydrogen 
content within the limits of accuracy of the 
parameter determinations. The lines in the 
pattern for vacuum-distilled manganese are quite 
sharp, and the K-alpha doublets are resolved for 
diffraction angles greater than 45°. The parameter 
could thus be determined to a greater accuracy 
than for the other samples, where the lines were 
broader and none of the doublets were resolved. 


EFFECT OF DEGASSING TREATMENTS ON 
LATTICE PARAMETERS AND 
DIFFRACTION PATTERNS 


Several batches of electrolytic manganese were 
degassed, according to the procedure suggested in 
a previous publication,® to determine the effect of 
degassing om the parameter. The metal was 
heated to temperatures from 300° to 925°C at 
reduced pressure until the pressure again de- 
creased to its initial value before heating.The 
treatment times were about 1 hour for tempera- 


TABLE II. Average parameter for various 
degassing treatments. 














Material Parameter 
Original manganese metal 8.904 
Manganese degassed at 300°C 8.902 
Manganese degassed at 500°C 8.8944 
Manganese degassed at 700°C 8.8948 
Manganese degassed at 925°C 8.8935 
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tures up to 700°C and one-half hour at 925°C. 
This did not remove all the gas, so samples of 
these batches were tested for gas content by 
heating them to 500°C in a hydrogen atmosphere 
at normal atmospheric pressure and determining 
the amount of gas liberated or absorbed by the 
sample in reaching equilibrium under these con- 
ditions. The solubility of hydrogen in manganese 
under these conditions is known and the amount 
of gas that was in the metal originally was either 
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Fic. 1. Variation in average parameter with hydrogen gas 
content of the metal. 


greater or less than this value by the volume 
liberated or absorbed. All the samples showed 
less than 10 cc of gas per 100 grams of metal or 
less than 4 percent of the original amount. 
X-ray diffraction patterns of each sample were 
then made, and the parameters were determined. 
The results are shown in Table II as compared 
with the parameter for the original metal. 

Removing the gas from the metal has little, if 
any, effect on the parameter at temperatures up 
to 300°C. At 500°C and over, however, the 
parameters return nearly to hormal value. The 
parameters listed for these three samples are 
essentially equal (within experimental error, 
which is about 0.0005A), and the average value, 
8.8942A, agrees well with the published value. 

A better comparison can be made from the 
diffraction pattern shown in Fig. 2; Figs. 2a, 2b, 
2c, 2d, and 2e show the patterns for the samples 
treated at 925°, 700°, 500°, 300°C and the normal 
untreated metal, respectively. The patterns on 
the 500°, 700°, and 925°C samples all show sharp 
lines with the K-alpha doublets resolved. The 
925°C pattern seems to be the sharpest; the 
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a. b. c. 
Metal degassed Metal degassed Metal degassed 
at 925°C at 700°C at 500°C 


g. 

Electrolytic 

i. manganese 

e. Metal trans- containing 
Normal formed from 615 cc of He 


d. 
Metal Suge electrolytic gamma to per 100 g of 
at 


manganese alpha manganese 





Fic. 2. X-ray diffraction patterns. 


500°C pattern is definitely broader. Densitometer 
measurements shown in Table III show the 925° 
and 700°C patterns to be essentially the same, 
though visually the 700°C pattern would seem to 
be somewhat broader. The 300°C pattern cannot 


be distinguished from that for the untreated 
metal, showing that the metal is essentially the 
same even though the gas has been removed. 
The diffraction patterns also indicate other 
variations in the structure aside from the differ- 
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ences in parameters of the various samples. The 
patterns shown in Figs. 2a, 2b, and 2c show sharp 
lines, while those in Figs. 2d and 2e show broad 
lines. Figure 2f shows the pattern for manganese, 
which was plated as flexible gamma-manganese 
and transformed to alpha-manganese, while 
Fig. 2g shows the pattern for a sample of manga- 
nese containing an unusually large amount of gas. 
Both of these patterns also have broad lines. 

Two factors cause lines in a diffraction pattern 
to be broad: (1) variations in the parameters of 
the crystals constituting the individual sample or 
(2) variations in the size of the crystal aggregates 
in the different samples. Parameter variations in 
the crystals constituting the sample would obvi- 
ously cause broad lines because the resulting 
pattern would be a composite of patterns from 
the individual crystals, and the lines in the indi- 
vidual patterns would be displaced relative to 
each other by amounts depending on the varia- 
tions in the parameter. Measurements made at 
the densest portion of the lines determine the 
average parameter, and the width of the lines is a 
measure of the deviation from the average in the 
respective crystal units in the specimen. 

Variations in the size of the crystal aggregates, 
which act as units to the x-ray beam, influence 
the width of the diffracted beam and the width of 
the diffraction lines. This phenomenon has been 
made the basis of a method for estimating the 
size of the crystal aggregates which act as units 
to the x-ray beam.’ In the present case both of 
these factors may be contributing to the observed 
line width, so an attempt was made to determine 
the actual deviations in parameter. In the follow- 
ing discussion ‘‘deviation in parameter” will be 
used to mean the difference between the parame- 
ter of an individual crystal and the average 
parameter for a sample while “variation in 
parameter”’ will mean the difference between the 
average parameter of a sample and the normal 
parameter for manganese metal. 

X-ray data alone are not sufficient to separate 
the two effects, so the specimens were examined 


7 (a) R. Brill, “X-Ray Determination of Size and Shape 
of Submicroscopic X'tals,"" Kolloid Zeits. 55, 164-169 
(1931); (b) G. H. Cameron, and A. L. Patterson, ‘‘The 
X-Ray Determination of Particle Size,"" Symposium on 
Radiogfaphy and X-Ray Diffraction, A.S.M., p. 324 
(1936); (c) H. Mark, ‘“‘The Determination of Particle Size 
by the Use of X-Rays,"’ Trans. Faraday Soc. 25 (1929). 
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microscopically to detect any variation in crystal 
or grain size in the different samples. No differ- 
ence could be seen between the regular elec- 
trolytic manganese metal and the degassed metal 
for temperatures up to 700°C even though the 
diffraction patterns showed marked changes. In 
each case the grains were very small. At 925°C 
some evidence of larger grains was seen, but the 
diffraction patterns show little difference between 
the 925° and the 700° samples. It would seem, 
therefore, that crystal aggregate or grain size has 
little influence on the line widths and that the 
grains are all sufficiently large to produce sharp 
diffraction lines. We can then attribute the 
variations in the line widths to deviations of 
parameter only and estimate the magnitude of 
this deviation. 

The width of the lines in six of the patterns in 
Fig. 2, measured with a densitometer, are shown 
in Table III. 

The line widths were measured at an intensity 
equal to one-half the maximum. The sharpest 
lines were obtained on the samples degassed at 
700° and 925°C, the widths being essentially the 
same. The pattern for vacuum-distilled metal is 
not shown because it was identical with that for 
the metal degassed at 900°C. These three ma- 
terials should be the most uniform, in structure 
and, since they gave almost identical patterns 
and line widths, it was assumed .that they most 
nearly represented material with a uniform 
parameter. The observed line width in the pat- 


_ terns for these samples is then due to the optical 


characteristics of the x-ray camera, and the size 
of the crystal aggregates in the specimens. These 
factors were assumed to be constant, and this 
minimum line width was used to correct the other 


TaBLE III. Line widths and parameter deviations in 
various samples. 








Line width Parameter 
Material in mm deviation 
Regular electrolytic manganese 0.84 +0.21 percent 
Electrolytic manganese degassed 
At 300°C 0.84 +0.21 percent 
At 700°C 0.52 0 
At 925°C 0.55 0 


Electrolytic manganese containing 0.61 +0.06 percent 
615 cc per 100 grams of Mn 

Electrolytic manganese trans- 0.99 +0.31 percent 
formed from gamma-Mn to 
alpha-Mn 
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values, the increase in width relative to this value 
being attributed to non-uniformity in the 
parameters. 

The maximum parameter deviation required to 
explain the broadening of the lines is only 0.3 
percent but, in this case, we are not justified in 
assuming that the crystal aggregates are the 
same size as in the other samples because the 
plating conditions required for producing gamma- 
manganese are different from those used for 
plating alpha-manganese. In the other samples a 
deviation of 0.2 percent in parameter value is 
indicated. It might be expected that the devia- 
tions would be greater with the larger amounts of 
gas, but this does not seem to be the case, since 
the samples containing the most gas show a 
deviation of only 0.06 percent in parameter. This 
observation leads to some possible conclusions 
regarding the solubility of hydrogen in manganese 
which will be discussed later. 


DISCUSSION OF RESULTS 


a. Variation of Average Parameter with 
Hydrogen Content 


The data in Table I and Fig. 1 indicate that the 
lattices of the manganese crystals are distorted 
by the hydrogen present in them. This distortion 
produces an enlargement of the average lattice 
parameter roughly proportional to the amount of 
gas in the metal or about 0.0003 percent per cc 
of gas per 100 grams of metal. The lattices are not 
necessarily distorted uniformly; any deviations 
from the average parameter should cause a 
widening of the lines in the x-ray diffraction 
pattern, the widening being proportional to the 
amount of the deviation. The lack of resolution 
of the K-alpha doublets in the patterns obtained 
from the samples containing gas indicates that 
such deviations may exist within the lattices of 
those samples, but the vacuum-distilled metal 
seems to have rather uniform lattice dimensions. 


b. Effect of Degassing Treatments on 
Lattice Parameters 


The lattice parameters obtained for the 
degassed manganese metal (shown in Table I!) 
and the x-ray diffraction patterns for these 
samples (shown in Fig. 2), show definitely that it 
is only necessary to remove the gas from the 
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metal and heat it to 500°C or higher to produce a 
metal with a normal and uniform lattice parame- 
ter. Since this treatment cannot affect any of the 
impurities, other than the amount of hydrogen 
present in the metal, the enlargement of the 
lattice parameter in the normal electrolytic metal 
must be due to the presence of the hydrogen 
introduced into the metal as atomic hydrogen 
during the plating process. The amount of hydro- 
gen is about 0.02 percent by weight; it occupies 
a volume at normal conditions twenty‘ times that 
of the metal containing it. It is not unreasonable 
to expect that this highly compressed gas would 
have some effect on the structure of the metal 
containing it. Degassing the metal at 300°C has 
virtually no effect on its lattices, probably be- 
cause the mobility of the atoms at this tempera- 
ture is too low to permit any appreciable read- 
justment of the lattice in the time required for 
removing the gas. 


c. Significance of Broadening of Lines in the 
Diffraction Patterns 


The line widths shown in Table III can be 
taken as reliable indices of the relative broadness 
of the lines in the various diffraction patterns, but 
the calculated deviations in parameter are sub- 
ject to too much error to be considered as more 
than rough estimates. It is evident, however, that 
the broadening of the lines observed in the 
diffraction patterns might be explained by a 
slight non-uniformity in the lattice parameters 
and that the non-uniformity is not proportional 
to the amount of gas in the metal. 

The average parameters for the regular elec- 
trolytic manganese, the manganese transformed 
from the gamma- to the alpha-form, and the man- 
ganese with the high hydrogen content were 
enlarged 0.11, 0.045, and 0.18 percent, re- 
spectively (see Table I). The deviations in 
parameter shown in Table III are larger than this 
and indicate that some of the crystals in the 
metal had parameters less than normal. This may 
be true, in which case the average enlargement 
must exceed the average compression of the 
lattice so that a net enlargement in parameter is 
observed. On the other hand, the calculated 
deviations may be too large because of improper 
measurement of the line width or improper cor- 
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rection for the width due to grain or crystal 
aggregate size. Possible errors from the first 
source cannot be discussed here because no ade- 
quate theory is available for determining the 
relation between grain size, parameter deviation, 
and line width. Errors from the second source can 
produce the effects noted if we assume that there 
is some grain growth in the manganese during the 
degassing treatments. The grains would be larger 
in the samples treated at the higher tempera- 
tures, and the lines in their diffraction patterns 
would be sharper. The correction for grain size, 
based on the line width of the 700° and 925°C 
samples, would be too small and we would be 
attributing too much of the widening effect in the 
other samples to parameter deviations. A more 
nearly correct value for the parameter deviation 
might then be obtained by considering that the 
deviation’ above and below the average parameter 
value is equal to the enlargement of the average 
parameter relative to that for normal alpha- 
manganese metal. 

The parameter deviations shown in Table III 
do not vary with gas content in the manner that 
might be expected. The deviation is zero at zero 
gas content, high at 95 cc per 100 grams of metal, 
and decreases with increasing gas content up to 
the observed maximum gas content of 615 cc per 
100 grams of metal. This behavior indicates that 
the deviation might again be zero at some value 
of gas content near 615 cc per 100 grams of 
manganese. In this case each crystal would be 
enlarged the same amount and would contain the 
same amount of gas. This would not be likely to 
occur unless the metal were saturated with gas 
and, in this case, there should be an integral 
number of atoms of gas in each unit cell of alpha- 
manganese. The number of unit cells of alpha- 
manganese in 1 gram of metal is proportional to 


ON X-RAY PARAMETER 29 


the reciprocal of the product of the atomic 
weight, the number of atoms in a unit cell, and is 
some constant times 0.3110-*. On the same 
basis, 615+50 cc of gas in the metal would 
provide (0.55+.05) X10-* atom of hydrogen, or 
nearly two atoms per unit cell of alpha-manga- 
nese. If we assume that each unit cell absorbs 2 
atoms of hydrogen, then it would require 695 cc 
of hydrogen per 100 grams of manganese to 
saturate the metal, and only about 90 percent of 
the unit cells in the sample containing 6.5 cc of 
gas per 100 grams of metal would be saturated. 
The remaining 10 percent of the unit cells would 
be distorted various amounts and cause the 
diffraction lines to be wider than if the sample 
were uniform. The deviation shown in Table III 
is consistent with this conclusion, and it would be 
anticipated that slightly higher gas contents 
would produce sharper diffraction lines. 

No samples of manganese were found contain- 
ing more than 615 cc of hydrogen, so we have no 
experimental evidence to show whether or not 
manganese will absorb more than 2 atoms of 
hydrogen. The fact that the regular electrolytic 
metal, containing 250 cc of gas per 100 grams, 
shows definitely fuzzy lines, indicates that the 
unit cells take 2 rather than 1 atom, however, 
because if the latter were true, the unit cells 
could be uniform with about 350 cc of gas per 100 
grams of metal, and the observed diffraction 
patterns show no evidence to support this 
conclusion. 
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Electron temperatures in the positive column and in the Faraday dark space of a cold cathode 
glow discharge have been measured by the probe method when the gas pressures of dry air and 
carbon monoxide were varied in the range from 0.06 to 1.2 mm Hg. For the fast group, electron 
temperatures decrease with pressure while for the slow group they are independent of it. The 
electron current densities increase and the space potentials decrease when the gas pressure is 
increased. Different current-voltage characteristic curves have been found in the positive 
column and in the transition region from the negative glow to the Faraday dark space. Theo- 
retical interpretations of these different characteristics could be obtained either from Langmuir’s 
theory of probe or from Druyvesteyn’s theory of electron energy distribution by assuming the 
existence in the boundary between negative glow and Faraday dark space of a drift fast electron 
beam superimposed upon an isotropic slow electron distribution. Spectroscopic investigations 
of the bands of the neutral and ionic molecules of Nz and CO in different sections of discharges 
also support the supposition. 



























1, INTRODUCTION Probe measurements'*" indicate that fast 
electrons with energies of the order of that ob- 
tained by passing through the cathode fall space 
are present in the Faraday dark space and some- 
times in the positive column and that they are 
distributed according to the Maxwell-Boltzmann 
law. Some authors,'® ® however, suspect that the 
primary electrons from the cathode dark space 
are able to penetrate far into the negative glow 
or even into the Faraday dark space, and it was 
suggested that these electrons are formed by 
some mechanism of the converse of ionization by 
collision, when a system of two slow electrons and 
one positive ion pass into a unit consisting of an 
excited molecule and a single fast electron. In a 
neon glow discharge Druyvesteyn™ has observed 
spark lines with high excitation energies, which 
proves that electrons exist with an energy corre- 
sponding to more than 80 percent of the cathode 
fall. But the concentration of these electrons 
could not be determined. A direct determination 
of the energy of the electrons'*-® showed that for 
a cathode fall above 400 volts, most electrons in 





UR knowledge of the mechanism of the glow 

discharge has not been advanced greatly, 
partly because the fundamental processes are 
only in a small part quantitatively known, and 
partly because of mathematical difficulties!“ 
which arise in the description of the observed 
phenomena of the discharge. But it has been 
established that, in order to maintain a self- 
sustained electrical discharge, the electrons gener- 
ated in the cathode region are mainly due to the 
y-process, which gives electron liberation at the 
cathode by positive ion impact.5-* These elec- 
trons are accelerated in the cathode fall space and 
move with nearly a homogeneous velocity 
towards the anode.* Of course, most of the 
primary electrons are scattered in their paths and 
retarded by collisions with gas molecules. The 
light emission in the negative glow is probably 
due to the result of the excitation of the gas 
molecules by these fast electrons entering this 
section of discharge from the Crookes dark space. 




















1P. M. Morse, Phys. Rev. 31, 1003 (1928). amas. 

2 J. J. Thomson, Phil. Mag. 8, 393 (1929). 10K. G. Emeleus and W. L. Brown, Phil. Mag. 7, 17 
3 W. Rogowski, Archiv. f. Elektrotechnik 26, 643 (1932). (1929). 

*W. Weizel, R. Rompe, and M. Schon, Zeits. f. Physik 1K. G. Emeleus, W. L. Brown, and H. M. Cowan, Phil. 












112, 339 (1939); 113, 87 and 730 (1939). Mag. 17, 146 (1934). 
5H. A. Wilson, Phys. Rev. 8, 227 (1916). 122K, G. Emeleus and O. S. Duffendack, Phys. Rev. 47, 
*A. K. Brewer and R. R. Miller, Phys. Rev. 42, 786 460 (1934). 
(1932). 13M. J. Druyvesteyn, Zeits. f. Physik 62, 764 (1930); 
7L. B. Loeb, Rev. Mod. Phys. 8, 267 (1936). Physica, 1, 427 (1934). 
8’ M. J. Druyvesteyn and F. M. Penning, Rev. Mod. 144A. K. Brewer and J. W. Westhaver, J. App. Phys. 8, 
Phys. 12, 87 (1940) ; also references given there. | 779 (1937). 
* J. J. Thomson, Phil. Mag. 48, 4 (1924). 18]. F. Lehmann, Proc. Roy. Soc. A115, 624 (1927). 
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ELECTRON TEMPERA 





Fic. 1. Variation of tube current 
with gas pressure. 














the negative glow had an energy which was 
almost equal to the cathode fall. Determinations 
of this kind, however, are not very reliable as a 
third electrode which disturbs the discharge must 
be brought into the negative glow. In the present 
investigation it is attempted to show the exist- 
ence, in the transition region from negative glow 
to the Faraday dark space, of primary electrons 
coming from the cathode dark space. Probe 
measurements have been made in the positive 
column and in the Faraday dark space when the 
gas pressures in the discharges were varied. 
Different current-voltage characteristic curves in 
these two regions have been found. Theoretical 
interpretations of these different characteristics 
could be obtained either from Langmuir’s 
theory'*"? of probe or from Druyvesteyn’s 
theory'® of electron energy distribution by as- 
suming the existence in the boundary between 
negative glow and Faraday dark space of a drift 
fast electron stream superimposed upon an 
isotropic slow electron distribution. Spectroscopic 
investigations of the bands of the neutral and the 
ionic molecules of N2 and CO in different sections 
of discharges also support the supposition. 


161. Langmuir, Phys. Rev. 26, 585 (1925). 

‘7H. M. Mott-Smith and I. Langmuir, Phys. Rev. 28, 
727 (1926). 

1®M. J. Druyvesteyn, Zeits. f. Physik 64, 781 (1930); 
Physica 10, 69 (1930). 


2. APPARATUS 


The form of the discharge tube used was 
cylindrical, 2 cm in diameter and 10 cm in length 
with two aluminum disk electrodes of 1.2 cm 
diameter sealed through a CR lead-in wire 
coaxially on each end of the tube at a separation 
of about 7 cm. The probe was a 6-mil tungsten 
wire 6 mm long and spot-welded to a lead-in wire 
mounted in a small glass tube fused perpen- 
dicularly to the discharge tube at a distance of 
about 1.5 cm from one of the disk electrodes. The 
discharge tube was sealed by means of a side- 
tube on to the vacuum system made of soft glass 
tubings, and pumped by a Cenco steel mercury 
diffusion pump backed by an oil pump. The 
pressures in the tube were measured by two 
McLeod gauges for different ranges. The dis- 
charge tube and electrodes were degassed by 
running a discharge with an induction coil and 
pumping out alternatively over days. 

The gases used were dry air and carbon mon- 
oxide. The latter was generated by heating oxalic 
acid in the presence of concentrated sulfuric acid, 
carbon dioxide gas and water vapor being re- 
moved by lime water and anhydrous phosphorus 


pentoxide, respectively. Before being employed 
to take readings the vacuum system and tube 
were flushed a few times with the gases. Because 
the gas pressure was not too low and the dis- 
charge temperature was not more than 35°C, it 
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Fic. 2. Variation of electron 
temperature with pressure. Upper 











curves, fast group; bottom curves, 
slow group. 














Fic. 3. Variation of electron current 
density with pressure. 








was thought not worth while to use liquid air to 
remove any residual mercury vapor in the tube 
although traps were provided on the vacuum 
system. 

The discharge was excited by a 500-volt 
rectifier stabilized by neon bulbs, with one of the 
electrodes grounded. The wiring in the probe 
circuit was arranged in the usual way. The 
potential of the probe was varied and the current- 
voltage characteristic values were obtained by 
means of a microammeter in combination with a 
universal shunt so that the current could be 
measured conveniently from a few microamperes 
to a few milliamperes. 


02 wm Mig 


3. EXPERIMENTAL RESULTS 

In the gas pressure range between 0.06 and 1.2 
mm Hg with a potential of 500 volts applied 
across the discharge tube, the tube current was a 
few milliamperes and so the discharge was 
probably of the normal glow type. On the anode 
side two or three striations could usually be 
observed. They were particularly sharp in the 
case of carbon monoxide gas. As the gas pressure 
was increased in the above-mentioned range, the 
discharge current increased linearly in the case of 
CO, but for air after 0.5 mm Hg the current 
seemed to approach a constant value as shown 


in Fig. 1. 
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Fic. 5. Variation of electron concen- 
tration with pressure. 











When the probe was at a very negative po- 
tential with respect to the anode, only a positive 
ion current was recorded. By decreasing the 
retarding potential the positive ion current at 
first diminished linearly and then rather rapidly, 
because some electrons of higher velocities were 
able to penetrate through the ion sheath around 
the probe and then were collected by the latter. 
As the retarding field was made still weaker, 
electrons of lower velocities could also be col- 
lected by the probe. Then an electron current 
became appreciable and continuously increased 
until one came to a point where the potential of 
the probe was the same as that of the surrounding 


1 
CE tn. te Hg 


ionized gas. In analyzing the characteristic data, 
a large scale plot of positive ion current versus 
probe potential was made. The extrapolation of 
the straight line tangent to the ion current to the 
smaller negative voltages gave the contribution 
of the positive ions to the current, and the 
electron current could then be obtained from the 
difference between this and the total observed 
current at that potential. In the regions of the 
positive column and the Faraday dark space near 
the anode side in the discharge, typical semi- 
logarithmic plots of electron current against 
retarding potentials are similar to the B curves in 
Figs. 6 and 7. 
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Fic. 6. Current-voltage char- 
acteristics in air discharge at a 
gas pressure 1.25 mm Hg, tube 
voltage 490 volts, tube current 
4.5 ma. Curve A in Faraday 
dark space, curve B in positive 
column. 
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If the velocities of the electrons were described 

by a Maxwell-Boltzmann distribution, the semi- 
logarithmic plot would be, according to the probe 
theory of Langmuir and Mott-Smith,’’ a straight 
line, of which the slope determined a temperature 
corresponding to the electron velocities. At higher 
retarding potentials the observed deviation from 
linearity indicates that the distribution of electron 
velocities in the discharge must be rather compli- 
cated. However, one might assume that there 
existed more than one group of electrons, each 
of which was distributed according to Maxwell's 
law. The slope of the straight line drawn along 
the higher retarding potential characteristic 
curve gave the temperature of the fast electron 
group, and the slope of another straight line 
determined from the differences between the 
extrapolation of the first straight line and the 
observed values gave the temperature of the slow 
electron group. Sometimes a third electron group 
was observed. In the case of a discharge in CO, 
the second group came out quite definitely and 
was analyzed. In the case of air, the presence of 
the second group was revealed only at higher 
pressures. 


197, Langmuir and H. M. Mott-Smith, Gen. Elec. Rev. 
27, 449, 538, 616, 762, and 810 (1924). 











Figure 2 gives the data obtained for the varia- 
tion of the electron temperature with pressure. 
For the fast group, electron temperatures de- 
crease strongly when the pressure was increased, 
while for the slow group, the electron tempera- 
ture is independent of the pressure. 

Figure 3 shows the data for the variation of the 
electron current density with pressure. 

Figure 4 presents data for the change of space 
potentials at a point in the discharge when the 
pressure was increased. 

Electron concentrations in the ionized gas can 
be calculated from the relation 


n= (I/e)(2am/kT)}, 













where e and m are the charge and mass of the 
electron respectively, k is the Boltzmann con- 
stant, J is the electron current density, and T is 
the electron temperature. The results of the 
calculation are plotted in Fig. 5. 

Figures 6 and 7 are the current-voltage charac- 
teristic curves in air and CO respectively for 
discharges in the Faraday dark space near the 
boundary with the negative glow (curves A) and 
in the positive column (curves B). These two 
regions yielded quite different characteristic 
curves. Near the space potential on the charac- 
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teristics we had a straight line portion in the 
Faraday dark space, which determined a temper- 
ature if the part of the characteristic at the higher 
retarding potentials was neglected. In the case of 
air at a pressure of 1.25 mm Hg this determined 
temperature is higher than the temperature of 
the fast group in the positive column, while in 
CO at a pressure of 0.87 mm Hg the determined 
temperature is lower than that in the positive 
column. 


4. INTERPRETATION OF THE VARIATION OF ELEC- 
TRON TEMPERATURES AND SPACE 
POTENTIALS WITH PRESSURE 

The observation of the decrease of electron 
temperature with pressure is quite similar to the 
measurements of Seelinger and Hirchert?® in the 
positive column of neon and argon discharges, of 
Groos” in argon, and of Klarfeld” in mercury. As 
the pressure increases, there are more chances for 
the electrons to impact with gas molecules and 
consequently the mean kinetic energy of electrons 
will be diminished because of energy loss. Hence 
the electron temperature decreases. The slow 
electrons are formed from primaries after being 
scattered and retarded by numerous collisions 
with gaseous molecules. Since these ultimate 
electrons have already suffered many elastic 
impacts, further change of pressure will not 
probably influence their mean kinetic energy. 
Thus the slow electron temperature keeps con- 
stant with pressure variation. 

By comparing the values given in Fig. 2 with 
that in Fig. 4 for any given gas pressure, one finds 
that the electron temperatures decrease with 
space potential, because in weaker fields electrons 
gain lower energy. This is in agreement with the 
observations of Duffendack and Chao* who 
found that electron temperature diminishes as 
one goes from cathode to anode in a nitrogen 
glow discharge. 
5. TRANSLATIONAL AND RANDOM 
ELECTRON VELOCITIES 


The difference of the volt-ampere characteristic 
curves obtained at the junction of the Faraday 


1931) Seelinger and R. Hirchert, Ann. d. Physik 11, 817 

*1 OQ. Groos, Zeits. f. Physik 88, 741 (1934). 

2 B. Klarfeld, Tech. Phys. USSR 4, 44 (1937); 5, 725 and 
919 (1938). 

%Q.S. Duffendack and K. T. Chao, Phys. Rev. 56, 176 
(1939). 
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dark space with the negative glow and in the 
positive column does not seem to be explained 
by merely dividing electrons into isotropic 
Maxwellian groups, because it is not reasonable 
to neglect the electron current for higher retard- 
ing potentials in the characteristic curve from the 
dark space. These characteristics, which are con- 
cave downwards for a considerable part in higher 
retarding potentials, are different from the 
anomalous curves discussed by Emeleus and 
Brown,” and must be due to some causes other 
than positive space charges. 

Let us assume that the resultant velocity 
distribution of the electrons near the Faraday 
dark space can be analyzed into a part with uni- 
form translational motion and a_ superposed 
Maxwellian part with temperature motion, and 
that the velocities due to isotropic temperature 
motion are small compared to those of transla- 
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Fic. 7. Current voltage characteristics in CO discharge 
at a gas pressure 0.87 mm Hg, tube voltage 520 volts, tube 
current 3.4 ma. Curve A in Faraday dark space, curve B in 
positive column. 


24K. G. Emeleus and W. L. Brown, Phil. Mag. 22, 898 
(1936). 
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Fic. 8. Inverse probability integral against probe potential. 


tion. The theory of the current to a cylindrical 
probe with its axis perpendicular to the direction 
of a stream of the drift electrons superimposed on 
a Maxwellian distribution has been treated by 
Mott-Smith and Langmuir,'’ and the result is 
expressed in a quite complicated series. However, 
in the region where the retarding potential of the 
probe is less than the equivalent voltage of the 
drift electrons, the volt-ampere characteristic of 
a cylindrical probe is nearly the same as that of a 
plane probe placed perpendicularly to the direc- 
tion of the drift velocities. Langmuir'®** has 
shown that if Jo is the current when all the 
electrons are collected by the probe, the electron 
current on the probe is given by 


T= (Io/2)[1+P(—A) ], (1) 


where P(—A) stands for the probability integral 
—x 
P(-2)=(2/v) f exp (—x?)dx. (2) 
0 


Here X is proportional to the difference of the 
velocity of the electron in the mass-motion and 
the velocity component of any considered elec- 
tron in the same direction, or defined by the 


relation 
A=(V Vi-V/. Vp) /V Vu. (3) 


Here V, is the retarding potential of the probe 
with respect to the ionized gas (or the sum of the 
potential of the probe with respect to the ionized 
gas and of the space potential), V; is the potential 


*6 1, Langmuir, Zeits. f. Physik 46, 271 (1928). 
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correspending to the drift velocity of the primary 
electrons, and Vy is the potential correspending 
to the temperature motion. If V; and V, are not 
too far from each other, the geometrical mean 
(V.V,)! can be replaced by the arithmetical 


mean 3(V;+V,). Equation (3) becomes 

(Vit Vp) 

A=(Vi/ vay [I], 
(2(V.Vm)!) 


Thus d is a linear function of the probe potential. 
From Eg. (1) and Eq. (2) we have 


P-'((2I —Ipo)/Io) =constant+ V,/2(V;V)!. 


(4) 


If the plot of the inverse probability integral 
against the probe potential is a straight line, the 
supposition of the existence of the superposition 
of the drift electrons upon a Maxwellian group is 
confirmed. The intersection of the straight line 
with the zero line determines the average energy 
of the drift electrons and the slope gives that of 
the random electrons. Figure 8 is such a plot 
from Eq. (4), and the results show that the ob- 
served values fall fairly well on straight lines and 
thus the experiments justify our separation of the 
electrons into translational and random parts. 
In the case of air the straight line intersects the 
zero line at 68 volts.”® This is the average energy 
of the electron beam, while the value directly 
determined from the curve A in Fig. 6 is only 50.7 
volts, since in the latter method of determination 
the higher energy electrons have been neglected. 
From the slope of the straight line we have 


2(V:Vu)'=44.8. (S) 


The cathode fall potential was not measured, and 
the fastest electrons that were collected by the 
probe have an energy about 102 volts, which 
could not be accurately determined. If this value 
is taken as the energy of the electron beam, from 
Eq. (5) we have Vy=4.93 volts, which de- 
termines the temperature assigned to the random 
electrons. In the positive column the slow elec- 
trons have a temperature corresponding to 3.34 
volts. Thus these two methods give about the 
same value of temperature for the slow electrons. 

In the case of CO, the determined average 

26 The value thus determined may not be very accurate 


because of use of the theory of current-voltage charac- 
teristics for a plane collector. 
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energy of the drift electrons is 65 volts and the 
directly measured value from curve A in Fig. 7 
is only 13.4 volts. The energy of the fastest elec- 
trons that reach the probe is 35 volts, while the 
average energy is 65 volts, so the primary 
electrons can not have energies less'than this. If 
Vi=65 volts and 2(V,V,)'=15.2, one finds 
Vw =0.89 volt. In the positive column the slow 
electrons have a temperature of 2.16 volts, which 
is a little higher. 


6. ENERGY DISTRIBUTION OF ELECTRONS IN 
POSITIVE COLUMN AND IN FARADAY 
DARK SPACE 


If it is assumed that an electron with a given 
velocity has angles of incidence upon the surface 
of the sheath of the probe less than a limiting 
value determined by the retarding potential, the 
energy distribution of the electrons can be 
derived. Druyvesteyn'® has obtained a distribu- 
tion function in the form 





=—(V,—V,)—., 


ee a= ” 4m di 
p ——_—— 
eA dV,? 


m 


where e and m are the charge and mass of the 
electron, A the surface area of the probe, 7 the 
electron current, and V, and V, are the gas and 
probe potentials respectively. Hence the measure- 
ment of the curvature of the current-voltage 
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Fic. 9. Observed electron energy distributions compared 
with Maxwellian distributions. Curve A for Faraday dark 
space, curve B for positive column. 
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Fic. 10. Observed electron energy distributions com- 
pared with Maxwellian distributions. Curve A for Faraday 
dark space, Curve B for positive column. 


characteristics will yield the distribution func- 
tion. The double differentiation from the charac- 
teristic curve, however, presents difficulties due 
to inevitable slight experimental errors. To 
simplify. the operation of the double differentia- 
tion, Sloane and Emeleus,”’ by making use of a 
mathematical identity, have obtained a modified 
formula 


2e(V,—Vp)} 
i. 
m 


4m, dlogi\? d*logi 
aoe Fy v>)| ( ) + | (6) 
e’A dV, dV,’ 





In a semi-logarithmic plot the current-voltage 
characteristic is straightened and the determina- 
tion of the second derivative is much facilitated. 
In the case where the semi-logarithmic plot has a 
linear portion, the second derivative of log 7 with 
respect to V, vanishes, while the first derivative 
is a constant. 

Curves B and A in Figs. 9 and 10 are the distri- 
bution functions of electrons calculated from 
Eq. (6) in the positive column and in the Faraday 
dark space respectively. In the positive column of 


27 R. H. Sloane and K, G, Emeleus, Phys. Rev. 44, 333 
(1933). 
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air discharge, the number of electrons observed 
with energy greater than 10 volts is less than that 
calculated from a Maxwellian distribution. Simi- 
larly in the CO discharge, for energies greater 
than 6 volts, the observed and calculated distri- 
bution functions begin to diverge. These observa- 
tions are similar to those of Druyvesteyn,'* and of 
Sloane and Emeleus”’ in the positive column of 
low voltage arcs in neon and argon. The observed 
number of electrons with about 14 volts in air or 
with about 11 volts in CO is particularly smaller 
probably because the energies of those electrons 
have been used up to excite the Nz or CO mole- 
cules, whose observed excitation potentials** have 
values mostly less than those voltages. 

In the Faraday dark space the distribution 
functions have a quite different form. For higher 
energies the number of electrons observed is so 
large that the Maxwellian distribution is negli- 
gible. These large values correspond to the points 
in the semi-logarithmic plot, where the portion of 
the curve has a curvature concave to the axis of 
voltage. The results of the distribution plot again 
show that in this section of discharge a fast 
electron stream is superimposed upon a small 
random electron distribution. The presence of 
this fast electron stream is of particular interest 
in the consideration of the excitation of molecular 
spectra in electrical discharges. 


7. ELECTRON ENERGIES AND 
MOLECULAR SPECTRA 


In order to get some information about the 


collisions of electrons with gas molecules, spectro- 


28 W. Jevons, Report on Band Spectra of Diatomic Mole- 
cules, pp. 78-82 (1932). 
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11. Spectrogram in air dis- 








charge. (a) in positive column, (b) in 
negative glow, (c) Fe arc. 





scopic studies were made on the regions of the 
negative glow near the Faraday dark space and 
the positive column in air and CO discharges. An 
image of the positive column in air was focused 
by means of a condensing lens on the slit of a 
Hilger E-1 spectrograph. The second positive 
bands of N2 are very intense on the photographs 
taken with Eastman Kodak plates at an exposure 
of about 70 minutes. When the negative glow was 
thrown close to the slit, the first negative bands of 
N.* are relatively stronger than the positive 
bands at an exposure of about 40 minutes, both 
bands being on the same plate in this spectral 
region as shown in Fig. 11. The excitation po- 
tential for bands of the second positive system is 
13 volts and the simultaneous ionization and 
excitation of the first negative system require 
19.6 volts.2® In the positive column the fast 
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gram in CO discharge 
(a) in positive column, 


(b) in negative glow, 
(c) Hg arc. 
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29 Recent studies of the excitation of Nz and N2* bands 
have been published in two papers by R. Bernard and R. C. 
Pankhurst, Science Abstracts 43, 917 (1941). These papers, 
however, are not available here. 
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electron group has an average energy of 16.6 
volts which is sufficient to excite the positive 
bands but not the negative ones. The negative 
bands can be easily excited in the negative glow 
by the fast electrons whose average energy 
measured is 68 volts. This seems to explain the 
different appearence of the spectra in these two 
regions. 

In case of CO a Hilger E-2 spectrograph 
happened to be employed to take pictures of the 
emission from the positive column ard the néga- 
tive glow. The spectrograms also show some 
difference in these two regions as seen in Fig. 12. 
The Angstrom bands of CO on both spectrograms 
are strong, but few of the comet-tail bands**-*! of 
CO*+ come to appear in the negative glow. 
Probably some of the Baldet-Johnson: bands are 
also present. The excitation potential of the 
Angstrom bands is 10.7 volts, while that for the 
comet-tail bands is 16.7 volts. In the positive 
column the fast electrons have an energy of 19.2 
volts, but in the region of transition from nega- 
tive glow to Faraday dark space they have 65 

*T. R. Merton and R. C. Johnson, Proc. Roy. Soc. 


A103, 383 (1923). 
* R. C. Johnson, Proc. Roy. Soc. A108, 355 (1925). 
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volts. It is obvious that the molecular ions are 
easily excited in the negative glow. That much 
higher energies are used to excite the bands of the 
molecular ions than that needed for the simul- 
taneous ionization and excitation, is because of 
the probability of excitation. Experimental in- 
vestigations of the electron-impact excitation 
function of the negative bands of N.* and the 
comet-tail bands of CO*+ have been reported by 
many people.*-* The conclusion is that the 
excitation functions for N2* and CO* molecules 
behave in a quite similar manner and both 
functions have a maximum at a potential about 
2.5 times the value of the excitation potential. 
The fact could be explained that though the 
excitation potential of CO* bands (16.7 volts) is 
less than that of N.* bands (19.6 volts), about the 
same amount of the electron energies is required 
to excite these bands. 

The authors wish to thank Mr. C. C. Yang 
who built a part of the vacuum system and 
Messrs L. Yang and F. S. Li of the Chemistry 
Department who helped to produce the carbon 
monoxide gas. 

% A. E. Lindh, Zeits. f. Physik 67, 67 (1931). 


33 R. Bernard, Comptes rendus 204, 488 (1937). 
%F. P. Bundy, Phys. Rev. 52, 698 (1937). 
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It has been shown that the mean number B(7) of discharges during any time T is given by 
B(T) = f"[1— pit) oeyat, 


where a(t) is the total density of radiation penetrating the tube at the time ¢ and p(t) is the 


solution of the integral equation: 


p(t) =f [1— p(x) Jote)d(e) 


for ¢>r and p(t) =1—e~* for OS tS 7. B(t) is explicitly obtained for the case that a(t) is con- 


stant. In this case it is also shown that 
B(T)—n 





aT= 


1—[B(T)—a}- 


|n| 


< (ar)? 


*1—(er)? ar<i. 


for 





INTRODUCTION 


HE errors in data obtained with G-M 

counter tubes, due to their finite resolving 

time have been subjected to mathematical analy- 
sis by several authors. 

In the past the expression m = mo exp [ —mor/T ] 
has been derived for correction of the number of 
discharges ‘‘n’’ obtained in time T with regard to 
the total number of particles mo penetrating a 
G-M counter tube and capable of producing 
discharges." 

For this derivation it was assumed that in 
order for a counter to register a particle, that 
particle should be preceded by a certain small 
time interval 7 during which no particle arrives 
at the counter. 

Most recently this simplified model was refined. 
It was assumed that a particle would produce a 
discharge at time (¢+dt) even if any number of 
particles arrived in the G-M counter tube during 
the time (t—r, ¢) provided, however, that these 
particles did not produce a discharge. 

The complete expression derived according to 
the second model for the number of particles 
missed is too complicated for practical use, 
especially if T is very large compared with 7, as it 
usually is. Consequently an asymptotic expres- 
sion for the evaluation of particles missed was 


1L. I. Schiff, Phys. Rev. 50, 88 (1936). 
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proposed by the same author.’? 
A(t) — nny? /(t+Nor). 


In addition the probability that m discharges 
will be produced in the G-M tube during the time 
T was derived on the basis that the probability 
of an electron passing through the tube in the 
time (¢, t+dt) is given by a(t)dt. The mean value 
and the variance of this probability distribution 
likewise were given. 

In the present paper it will be shown, with the 
second model, that the mean number of dis- 
charges, during the time 7, can be obtained from 
the solution of an integral equation. The special 
case when a(t) is constant is treated, and an ex- 
pression for the expected number of discharges is 
obtained for this case, which is somewhat simpler 
than the one previously derived.” 


COMPUTATION OF THE EXPECTED NUMBER 
OF DISCHARGES 

Let p(t) be the probability that a discharge 
will take place in the interval (¢—r,?¢). Let 
a(t)dt be the probability that an electron passes 
through the tube during the interval (¢, ¢+dzZ). 
The fact that p(t) is a continous function of ¢ is 
used here but its proof is postponed. In order that 
a discharge take place, it is necessary and 
sufficient that: 


2B. V. Gnedenko, J. Exper. Theor. Phys. U.S.S.R. 11, 
101 (1941). 
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1. An electron passes through the tube in the 
time interval (¢, t+dt). 

2. No discharge takes place in the time 
interval (t—r, ¢). 

Since 1 and 2 are independent events, the 
probability of one discharge in the time dt be- 
comes [1— p(t) Ja(t)dt. 

Then, since only one discharge can occur in the 
interval (t—r, ¢), the probability of one discharge 
in that interval is 


f [1—p(x)Ja(x)dx. ‘ 


But this is just p(t). 
Hence we obtain the integral equation 


p(t) = f [1 —p(x) Ja(x)dx. (1) 


Puttin 
~ + tape) glx), 


we then obtain 
t 
q(t) =1 -f q(x)a(x)dx. (2) 
For 0<t¢71; “i 


q(t) =exp |-f a(s)dx| 


Equation (2) may also be written in the form 
q (t) =qg(t—7)a(t—7) —q(t)a(t). (3) 


Let p(t) =1—g(t) be a solution of (1). The ex- 
pected number B(t) of discharges in the time ¢ is 
then given by 


B(t) -{ [1— p(x) Ja(x)dx. (4) 


0 


If t=sr+r’ with r’ <7, we may on account of 1 
write 


B(t)=¥ p(ar)+ f [1—p(x)Ja(x)dx. (5) 


a=1 


This reduces the problem to the solution of the 
integral Eq. (1). 
If a(t) =a independent of t, then Eq. (3) reads 


q (t) =aq(t—7) —agq(t). (3’) 


It is easy to see that Eq. (3’) together with the 
initial condition q(t) =e-*' for O< t <r determine 
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completely the solution of Eq. (3’). This solution 
is easily found as 


ima a'(t—ir)? 
q(t) = eatin 
i=0 1! 


for ar<t<(a+1)r. (6) 
Hence, if s is the largest integer smaller then 


t/r, 
ee EE 


i=0 4: 


e-**— #9) for ar Ct< (a+1)r, 


B(t)=> plar) +a f [1—p(x) dx 


a=1 


a=s i=a—l q iy (a —1) i 


=s5— > = e7e7(a—i) 
a=1 i=0 1! 
t ims a'(x—ir)! 
+af >. ————¢ «2 indy 
sr t=O 1! 


i=—s—1 qir* Bas—i 


ag~ 2 ~~ 2, em 


i-o0 4! p=1 


t ims a ‘(x —ir)é 
+af > ——e*- idx, (7) 


sr *=0 1! 





The number A(t) of particles which fail to 
produce discharges has previously been com- 
puted.? Clearly A (¢t)+B(t) =at. It was found that 


k=s a*(t—kr)* 





A(t) =ar >) ————e-t-#) 
k=1 (k—1)! 
k=s tr [t&—(k—-1)7 } 
k+1 lhe hes ie: SM 
oe) 
(k—1)! 


Xexp (—altk—(k—1)7 ])dty 


t (t,—sr)**! 
+are+ f Poaeer exp [—a(t,—sr) ]dty. 


8T Ss: 


In using the value B(t) which has been ob- 
tained in the present paper and the value for 
A(t) above, one obtains as a side result an 
interesting identity from the relation B(t) 
+A(t) =at. 
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To justify the above solution, it is necessary to 
show that p(t) is a continuous function of t. Now 
p(t+At)—p(t) is certainly smaller in absolute 
value than the probability of a discharge taking 
place during (¢, +At) plus the probability of a 
discharge taking place in the time (t—7, t—7r+ At). 
The probability of a discharge taking place 
during any time is smaller than the probability of 
an electron passing through the tube during this 
time. Hence 


t—r+At tat. 
| p(t+-At) —pt| <f a(s)de+ f a(x)dx. 


t—r 


Hence p(t) is continuous if 


t+At 
f a(x)dx 
t 
exists for all ¢. 


Formula 7 is hard to evaluate if S is large, as is 
usually the case in practical applications. For this 
case the asymptotic formula 


ny=n/(1—n7/t), (8) 


has been derived, where mp is the number of 
electrons passing through the counter and n the 
number of discharges during the time ¢. In this 
derivation it was assumed without proof that 
p(t) converges to a constant.” 

In the following a proof will be given for 
Eq. (8), for the case ar <1, and also an estimate 
for the error committed in using Eq. (8) for 
finite ¢. In the following, it will therefore be 
assumed that ar<1. The inequality ar <1 holds 
true in almost all practical applications, since r 
is of the order 10~‘, and ar is the mean value of 
the number of particles in the time interval r.* 4 


Putting 
ar 





p(t)= +e(t), 


ar 


then on account of (1) we have 


for t>r, 


e(t) = -af e(x)dx 


(9) 
for O<t¢€7; 





e(t) =1—e-** — 
1+ar 
>C.G. Montgomery and D. D. Montgomery, Phys. Rev. 


57, 1030 (1940); J. Frank. Inst. 231, 447, 509 (1941). 
4W. E. Ramsey, Phys. Rev. 57, 1022 (1940). 








AND H. B. MANN 
and 
ar 
le(t)| < for O<t<r, 
1+ar (10) 
|e(t)| <1 for all ¢. 


Applying the mean value theorem repeatedly, 
we have from Eq. (9) 


e(t) = —are(t— Or) =--- 
=(—ar)*e(t—aé’r), OC 6<1, 
< 
for ” 
argt¢(a+1)r. 


If t—aé’r is not yet in the first 7 interval, the 
process may be continued. If ¢(¢) is different from 
0, we must on account of | e(¢) <1, ar<1 finally 
obtain an equation of the form 


e(t)=(—ar)*e(t’) with O<t’ <1, ha. 
Hence on account of the relations (10) 
(ar)**! 


|e(t)| <——— for ar<t<(at+1)r. (12) 
1+ ar 


For the integral in Eq. (7) we have 





é a(t—sr) 
af (1-p)i=———"+5 (13) 
: as 1+ar 
where 
(ar)**? 
|5| < 
1+ar 


From Eqs. (7), (12), and (13) we have 


at 
B(t)=——+n, 
l+ar 
(14) 
(ar)? 
|n| <————[1 — (ar) **"]. 
— (ar)? 


The term 7 is very small even if ar is close to 1. 


Thus, for instance, if ar =0.95 from (14), |7| <10. 
Considering even that 7 is in the neighborhood of 
10-* then a should be in the neighborhood of 10°, 
and 7» is quite negligible indeed.* 4 

If the intensity of radiation is increased indefi- 
nitely, then m must approach t/r. Hence, 7 may 
be found by increasing the intensity of radiation 
until m remains constant. If 7% is this limit value 
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This is for all practical purposes equivalent to 
Eq. (8) as long as @r is not very close to 1. 


of n, then 


r=t/n. (15) 


Since for all practical purposes B(t) =n, at =m ACKNOWLEDGMENT 


there is obtained from Eq. (14) The grant of the Research Foundation of the 
Ohio State University and the support given by 
The Ohio State University Development Fund 


are greatly appreciated by the authors. 


u—n 


o= -_ (16) 
1 (n— (4/0) 
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The theory of volume effect in photoelectric emission is developed. Formulae are derived 
for the rate of electron excitation and for the photoelectric yield. The calculated threshold 
frequencies for volume effect are 5.9110" and 5.6910" for sodium and potassium respec- 
tively. The estimated photoelectric yields for these metals are of the same order of magnitude 
as those calculated for surface effect and are comparable with those observed experimentally. 
The volume effect should not, therefore, be neglected except in the immediate neighborhood 
of the threshold frequency. Approxmate estimation indicates that light absorption of sodium 
and potassium in the visible and ultraviolet regions should be largely due to quantum excitation 
of electrons. Accurate calculation of electron excitation and absorption requires detailed 


knowledge of electron wave functions. 


INTRODUCTION 


HE theory of photoelectric emission from 

metals is usually based on Sommerfeld’s 
model of free electrons for metals. Since a free 
electron cannot absorb the whole energy of a 
light quantum (the kind of absorption for photo- 
electric emission), it has been suggested that 
light absorption takes place only at the surface 
of the metal where the electrons are subject to 
the surface potential field. The phenomenon is 
then a surface effect. 

Wentzel' takes into account the damping of 
the light wave and accredits it with enabling 
the free electrons in the metal to absorb the 
energies of the light quanta. In this picture light 
absorption takes place inside the metal and is a 
volume effect. It has, however, been shown by 
Bethe? that the order of magnitude of emission 
given by this picture is too small. 

1G. Wentzel, Probleme der Modernen Physik (Sommer- 


feld Festschrift), p. 79. 
*H. Bethe, Handbuch der Physik 24/2, p. 472. 


The free electron model is too crude. Actually 
the electrons in a solid are in a periodic potential 
field. They can therefore absorb energies of light 
quanta independent of the damping of the light 
wave. Tamm and Schubin* made a rough esti- 
mation of the frequency, at which emission due 
to such absorption should begin, showing that 
it is much higher than the ordinary threshold 
frequency. Thereupon volume effect is neglected 
in the subsequent developments of the theory of 
photoelectric emission. It is our purpose to 
investigate this volume effect due to the periodic 
potential field inside the metal. 


DEVELOPMENT OF THE THEORY 


The potential energy of an electron in the 
crystal lattice of a solid is periodic 


Va), Vette, (1) 
The electron wave function is then of the form 


v= €*Tuz(r), 
+]. Tamm and S. Schubin, Zeits. f. Physik 68, 97 (1931). 
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where u,(r) is a function having the same 
periodicity as the potential. Optical phenomena 
with such electrons have been discussed in many 
places.‘ Let the vector potential of the light wave 
in the metal be 


nn a 
A =A; cos 2nr(t-"*) = £[ A pe? inalceeaivt 
HA pert inaleg 2rivt ), 


The damping can be neglected, since within the 
thickness in which we are interested, the wave 
is not appreciably damped. The perturbation 
due to this light wave is (eh/2rimc)A-VY¥. By 
using the method of variation of constants we 
get for an electron originally in the state k 











Qni 
ch exp =e Ex—hv)t—1 . 
Se bole Ev —E,—he Cw 
2ni 
exp (ey —E,.+hv)t-—1 - (2) 
* Ev —Ex+hy oy 


where 


Cue = J viaerinne. Vid; 


Com [vides Wad. 


As the wave-length of the light wave is much 
longer than that of the electron waves, we can 
neglect 2xnn/c compared with k and k’, then 


= ay 
Cvatie = af vi aly 
dg 
sad + IU, 
=a f —k) ry), alee Qz:Cx, x. (3) 


It has been pointed out by Kronig® that because 
of the periodic nature of the functions u(r), this 
integral is zero unless 


k’—k=2ng. 


This gives definite relations between the original 
state and the excited states; the electron cannot 
be excited from one state into any arbitrary 
state. 

*A. H. Wilson, The Theory of Metals, see references on 


p. 127. 
5 R. de L. Kronig, Proc. Roy. Soc. Al24, 409 (1929). 
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The probability of excitation of an electron 
in state k is 


S janacenl 
— | Dk, k+2ng | ~ 
> dt ‘ 
per unit time. The total rate of excitation, i.e., 
the total number of excited electrons per unit 
time, is 


Q d 
ie = ~ dt fff | Qk, k+20g | *dk.dk,dk, (4) 


where @ is the volume of the metal. The integra- 
tion is over all occupied states. At T=0°K the 
limits of integration are given by the surface of 
Fermi distribution. By changing coordinates, the 
integral can be written 


fff | Qk, k+2rg | *dudvdw, 


where u and v are coordinates in the surface of 
constant energy difference 


H= Exsorg—Ex—hv=0 (5) 


in the & space, and w is the coordinate perpen- 

dicular to this surface. Since d,4427, has a strong 

maximum at w=0 due to the factor 
exp (29i /h)(Exsong — FE, —hyv)t —1 


ie . ’ 


we may put the slowly varying factor 
Ck, k+209 = Ci, k+289 400 


in the integration with respect to w, and neglect 
the second term in d;,%427,. We then get 


eh \? 2 
[lau nsaeel*d0o= (<<) | @:|*| Cx, kong | w=) 
T 


x f /exp (24t/h)(Ex+2rq— Ex—hv)t— 1 4 


- Ww. 





The main contribution to this integral comes 
from near w=0. We may therefore extend the 
limits to — © and + and put 


Exsong—- Ex —hv= VHT po. (6) 


Carrying out the integration we get 


é 3 Ck k+2ng ia 
f [aan 2%0 |*dw = (=~) ht| Ge|? , 
2mc | VH | wad 





_ a, ak fe oe 
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Substituting this expression into (4) we get® 


Q 1 Cr,n+2n0|? 
yom ® (+ air f f 2d 
8x3 hoe ~ |vH| i) 


The integration is over the portion of surface 
H=0 bounded by the intersection of this surface 
with the surface of Fermi distribution. 

Not all excited electrons can be emitted. As 
the tangential components of k are continuous 
at the surface, on account of the continuity of 
wave function, the condition for an excited 
electron to be emitted is 


E(k. +2n., ky+2r7g,, k.+2ng:) 
2 


h 
= W+—_[(kyt+-2mgy)?+(ket+2ng.)*], (8) 
8x°m 


where W is the potential energy of an electron 
outside the metal. The volume of the metal Q 
being the product of the surface area and the 
thickness /, the emission current per unit surface 
area is 


J= 2eh ‘ee 3 @;|* 
2mc 





where the integration is over that part of the 
surface H7=0, which is bounded by the inter- 
sections of this surface with the surface of 
Fermi distribution and the surface (8). 


Effect of Surface Reflection 





in terms of the vector potential of the light wave 
inside the metal. To obtain the emission in 


terms of the incident light energy, we must find 
the relation between this vector potential and 
that of the incident wave. We can find this rela- 








2 sin 6; 
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The above equation gives the emission current * 
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tion from the classical electromagnetic theory 
by using optical constants of the metal: ” and k. 
The quantum theory of reflection developed by 
Schiff and Thomas’ shows that the classical 
relations hold true for the tangential com- 
ponents. For the normal components the relation 
differs from the classical theory only in the 
immediate neighborhood of the surface within a 
thickness of the order of 10-* cm. This is im- 
portant for the surface effect. For the volume 
effect we are considering it of no importance. 
Let the incident wave be plane polarized, and 
let the subscripts 7, r, and ¢ refer to the incident, 
reflected, and transmitted (refracted) waves re- 
spectively. 

(1) Electric vector of the incident wave 
parallel to the plane of incidence (x axis normal 
to the surface in the outward direction) : 


A;=(Qiz, — @;, cot 6;, 0) 


2aiv : 
Xexp ——(x cos 6;+y sin 6;+ct), 
c 


A,=(@,z, @rz cot 6;, 0) 


2riv : 
Xexp ——(—x cos 6;+y sin 6;+<ct), 
c 






A,= (Qéz, — Qiz cot 6:, 0) 








2rtv ned ; 
Xexp ——[(n—tk)(x cos 6:+-y sin 6;)+ct]. 
c 


The relation between @; and @; is given by 


(n—tk) sin 6,=sin 6,. 


Therefore 
; sin 6; 1 
sin 0,=———;_ cos 06, =———[ (n—ik)*—sin? 0; |. 
n—1tk n—1tk 


Now curl A=H and —(1/c)(dA/dt)=E. From 
the conditions of continuity for the tangential 
components H, and ¢,, we get 











a= 


Q; 
(n —ik)*+[(n—ik)*—sin? 0; ]*/cos 0; 


= @;Kuz, 
(10) 


2[ (n —ik)?—sin? 6; |} 


Qy= — @:, cot 6,.= —@ 


= @;Kuy. 





* This is substantially the same as Eq. (22.17), p. 457, 
Handbuch der Physik, 24/2. 





‘(n—ik)?+[(n—ik)?—sin? 6; ]!/cos 6; 


7™L. I. Schiff and L. H. Thomas, Phys. Rev. 47, 860 
(1935). 
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The transmitted wave in the metal is elliptically 
polarized since @:z and @y are complex with 
different arguments. For the special case of 
normal incidence 


Kuz=0; Kuy= —2/(1+-—ik). (11) 


(2) Electric vector of the incident wave per- 
pendicular to the plane of incidence: 


2riv ; 
A;=(0, 0, @i2) exp ——(x cos 6;+y sin 6;+ct), 
c 


2atv ; 
A,=(0, 0, @,z) exp ——(— x cos 0;+y sin 0;+ct), 
c 


wiv 





2 
A,=(0, 0, @:2) exp 


X[(n—ik)(x cos 6:+y sin 6,)+ct]. 


By the same procedure we get 
2 
‘1+[(n—ik)?—sin? 6,]!/cos 6; 
=@;:Ki. (12) 





Q:=@A:,.=@ 


For the case of normal incidence 
Ki=2/(1-+n—1ik). (13) 
The incident light energy per unit area of the 
surface per unit time is 
(axv?/2c) cos 6; | @;|?. (14) 
Substituting (10) or (12) into (9) and dividing 
by (14), we get the photoelectric yield, i.e., the 


charge emitted per unit area per unit incident 
energy 


_ Aceh (5. = Ef f 
2mc/ 8nx* cos 6; g 


| KuzC, nreeet KusCe, k+2rg | . 
| VH| 


4ceh “(5 ys 1 |K,|? 
= 2mc/ 8x* cos 6; 
| Ci, k+2ng|? 


x f ————dud?, 
~ | ViZ| 


for the two cases of electric vector parallel and 
perpendicular to the plane of incidence, respec- 
tively. Since the electrons excited deep inside 
the metal cannot be emitted because of the 








dudv, 
(15) 
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loss of energy by collisions with other electrons, 
1 in these equations is an effective thickness 
within which excited electrons reach the surface 
without much loss of energy. It depends, there- 
fore, on the collision process between electrons 
which is very difficult to calculate. Bethe’ 
estimates the effective thickness to be of the 


‘order of 10-7 cm. Rudberg’s experiments*® on 


inelastic scattering of electrons from _ solids 
shows that electrons of 100- to 200-volts energy 


‘can penetrate several atomic layers. Compton 


and Rose,® and Goldschmidt and Dember'® by 
experiments on photoelectric emission by light 
irradiating the back side of a thin layer, conclude 
that the electron mean free path in platinum is 
1.08 to 2.67 10-7 cm. It seems then reasonable 
to assume that / is of the order of 10-7 cm. 

To determine the photoelectric yield, we must 
know the electron energy levels for the calcu- 
lation of H and the electron wave functions for 
the calculation of Cy,%42*,. Such data have been 
obtained for a number of metals by the cellular 
method of Wigner and Seitz," and Slater." But 
the data have not been published in full. How- 
ever, for the determination of low frequency 
limits for light absorption and electron emission, 
we need to know only electron energy levels. 
Knowing the energy as a function of k, we shall 
have the equations of the three surfaces: (5), 
(8), and the surface of Fermi distribution with 
hv and g as parameters. The area, over which 
each of the integrals in the sums of (7) and (9) 
is to be integrated, can then be determined. For 
each integral (each value of g) there is a fre- 
quency below which the area and consequently 
the integral itself become zero. The frequency 
below which all integrals in (7) become zero is 
the long wave-length limit of light absorption; 
the corresponding frequency for (9) is the 
threshold frequency of emission. It has been 
shown by Wigner and Seitz, and Slater that 
except at energy discontinuities, the energy 


8 E. Rudberg, Phys. Rev. 50, 138 (1936). 

*K. T. Compton and L. W. Ross, Phys. Rev. 13, 295 
(1919). 

1H. Goldschmidt and H. Dember, Zeits. f. tech. 
Physik 7, 137 (1926). 

nu FE, Wigner and F. Seitz, Phys. Rev. 43, 804 (1933). 

2]. C. Slater, Phys. Rev. 45, 794 (1934). 

3 E. Wigner and F. Seitz, Phys. Rev. 46, 515 (Fig. 3), 
(1934). 

4 See reference 12 (Fig. 3). 
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Fic. 1 (a) Geometric representation of Eqs. (17) and (18). (b) Portion of the (u, v)-plane 


levels of sodium are given very closely by the 
expression for free electrons 


Ey = (h?/8x?m)k*. (16) 


We shall use this approximation for the deter- 
mination of threshold frequencies for sodium 
and potassium. 


THRESHOLD FREQUENCY 


Using the approximation (16) we find tha‘ 
Eq. (5) of the (u,v) surface becomes 


(h?/8x°m) (xg?+g-k) =hy, (17) 


which is the equation of a plane in k-space. The 
surface of Fermi distribution becomes a sphere 


k2+k/+ke =e. (18) 
The area over which we have to integrate for 
(7) is the portion of (17) bounded by its inter- 


section with (18). This is a circle with radius 
(Fig. 1a) 


[ , eet) 
r=| k,,?2— {| ——————_- 
g 





=(kn?—R2)?. (19) 


The area zr* is zero unless |R| <k» or 


2xrm 
TE — Bhim & Thy § wei + gk. (20) 


over which the integral of Eq. (9) is to be taken. 





This relation gives the low frequency and high 
frequency limits of absorption of light for the 
given g. The lowest of the low frequency limits 
for all possible values of g is the absorption limit 
of the metal. 

For electron emission the condition (8) must 
be satisfied, which by the approximation (16) 
becomes 


2 





(het 2ngs)* 2 W=E,+En 

2 
—(Re+km"), (21) 
“mM 


T 





where E, is the work function of the metal. This 
can be written 


kz > (ke?-+hm?)!— 2g, = ho. 


Let point A (Fig. 1a) be the origin of coordinates 
in the (uw, v)-plane and let the axis v be parallel 
to the (k,, k,)-plane. Then we have 


kz= jut+%Xo, 


where xo is the &,-coordinate of point A. From 
(17) we have 
xo = Rg. /g. 


j is cosine of the angle between axes u and k, 
j=[cos? (0A, k,)+cos* (OA, k.) }! 
=((gv/g)?+(gs/g)*}}. 







































Fic. 2. Values of a: 
(a) jho+ © (knt— ke)) <0 


(b) jho+™ (kn? — ket) >0 


(c) jho— (kn? — ket)*>0. 








The area { f{dudv over which the integrals in 
(9) are to be integrated is the part of (u, v)-plane 
(Eq. (17)) bounded by the circle wr? (see (19)) 
and the line 


u= (ko 





—Xo)/j 
shown by the shaded area in Fig. 1b. It is easily 
‘shown 


ff faudo=re ifa<c—l; 
ff feaudo=r {cos a—a(1 8), 


if —-l<a<+1; (22) 


ff faudo=o, ifad>+1; 


Ro— Xo _ ko— — Rgz/g 


le coaanemmens 
jr j(kn2— R28 

We see that a is a function of g and », R being a 
function of v. For each value of g we can deter- 
mine the limiting frequency », below which 
a>+1and { {dudv=0. The lowest of these for 
all possible values of g is then the threshold 
frequency. 

To determine », we shall solve the equation 


a=1, or (ko—Rg./g)?=7?(km?—R?). 
The two roots are 
R= (hogz/g) + j(Rm? —ko?)!. 


We note that g, must be positive since g.<0 
means decrease in the outward momentum, 
which certainly cannot give emission. There are 





and jko—" (knt— ke) <0 








then three possible cases 


(A) jko+ == (len?—ke2)! <0, 


g 






, gz , 
(B) JRo+—(Rm? — ko?) >0 
g 


and 


jo — (kn? — ke?) <0, 
g 


(C)  jko— (kn? — he?) >0. 
g 


It can be easily shown that the curves of a 
versus R for the three cases are as shown in 
Fig. 2. To have a<+1, R must observe the 
following conditions 


For case A: no limitations on R, 


For case B: R> = by— jk? —k?)!, 
g 


For case C: 
gz a . gz “a ¢ ° 
—ko— j(Rm? — ko”)! <R <—kot j(Rm? — ko?)!. 
g g 


Furthermore, if ko?>k»,? we have two possible 
cases 


D: ko<—km, 
E: ko>km.- 
We deduce from Fig. 2a and Fig. 2c: 


For case D: a always smaller than +1, there- 
fore no limitations on R. 


For case E: a always larger than +1, there- 
fore no emission. 


ahs a Ar ~_ - 





ko?) !. 


ssible 
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TABLE I. Low frequency limit of light absorption and the threshold frequency. 








Emission 1, 0,0 
surface 


Threshold 
frequency 


11, No, Na 


Threshold 
frequency 


"1, Ne, nN3 


Threshold 
frequency 


Ni, Ne, Ns 


12.6X 10" 
+1, 0,0 


5.91 10" 
+1, +1,0 


+1, +1,0 


+1, +1,0 


+1, 0, +1 


9.0 10" 


+i, —-1,0 


7.34X 10" 
0, +1, +1 


+1,0, +1 





Absorption limit 
Ni, Ne, Ns 


+1, +1,0 


3.14 10" 
+1,0, +1 


0, +1, +1 


Se 





Emission 1, 0,0 
surface 


Threshold 
frequency 


n1, ne, n3 


Threshold 
frequency 


Ni, Ne, N3 


9.57 X 10" 
+1, 0,0 


6.21 10" 
+1, +1,0 


+1, +1,0 


6.82 x 10" 


+1, —1,0 


+1, 0, +1 


Threshold 
frequency 


+1, +1,0 


Ni, Ne, Ns 


5.69 X 10" 


+1,0, +1 0, +1, +1 





Absorption limit 


1, Ne, Ns +1, +1,0 





4.96 x 10" 


+1,0, +1 0, +1, +1 








Case E gives a limit for g,, for it shows that a 
necessary condition for emission is ko <km 


(Rm? +k,?)'§—Rm 
2n 


This condition can also be directly deduced from 
(21). It gives the lower limit for g, and confirms 
our previous statement that g, must be positive. 
In addition to the limitations on R obtained for 
these different cases, we must also have | R| <km, 
otherwise r in (22) will be imaginary. It can be 
shown by combining the limitations on R that 
to have a<+1 we must have 





(km? +k,”)'—2xg, <Rm, £:> 


9 


h? 
hv >——(xg* — gkm) 
2rm 


if jko+ =< (kn?— ke?) <0, 
' (23) 


h2 
hvy> ——[ng*+geko oat ((g? —g,*) (Rm? - ko*))*] 
2xm 


if jko+ =< (kn? —ke?)'>0. 
g 


These relations give us the limiting frequency v,. 
The threshold frequency is the lowest of v, for 
all possible values of g. It can be shown that if 
g and g, were continuous variables, then the 
minimum hy, is 


(h?/8x°m)k,? = Ey, 


which is the ordinary expression for threshold 
frequency. The emission due to surface effect 
begins at this frequency. As g and g, can take 
only discrete values, the minimum », will be 
higher than this value. This confirms the general 
impression that emission due to volume effect 
begins at a higher threshold. In general », 
increases with g, but the lowest », does not 
necessarily correspond to the lowest value of g. 
Tamm and Schubin derived by simple reasoning 
an expression for v{,1/a),0,0) and assumed this to 
be the threshold frequency for volume effect. 
Actually, as shown later, for body-centered 
cubic lattices g does not have the value 
[(1/a), 0, 0]. But even if g has such a value, the 








50 as ws 








15 x10! 





Fic. 3a. Sodium: 


Curve 1. Experimental spectral distribution of photoelectric 
yield in coulomb/calorie. 

Curve 2. Ky X10’, plane (1, 0, 0). 

Curve 3. Ky X10’, plane (1, 1, 0). 

Curve 4. Ky X10’, plane (1, 1, 1). 


corresponding v, is not the lowest both for 
sodium and potassium as shown in Table I; the 
lowest v, corresponds to g=[(1/a), +(1/a), 0] 
and g=[(1/a), 0, +(1/a)] which have a larger 
value of g. 

The factor g comes from (1). For cubic 


lattices 
1 n 
g=-—(m, Ne, n3) a, 
a a 


where 1, m2, and m3 are integers, and a is the 
lattice constant. For body-centered and face- 
centered cubic lattices, the Fourier coefficients 
V, can be written: V,=A,S,, where S, is the 
structure factor. Sodium and potassium have 
body-centered lattices for which” 


S,=1+ cos r(m,-+2+n3). 


Thus S,=0 if (m;+m2+n3) is odd, and S,=2 if 
(m1+n2+n3) is even. Therefore for such lattices 
Viio0, Vary, ete., are zero, and g does not 
have the values [(1/a),0,0], [(1/a), (1/a), 


(1/a) ], etc. 


The maximum &k of the occupied states is 
given by 
km = (3x?)'(mo/v0)! 
where mo is the number of electrons per atom, 


and v is atomic volume. For metals with one 
conduction electron per atom (m)=1) and for 


% N.F. Mott and H. Jones, The Theory of the Properties 
of Metals and Alloys, p. 158. 
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Fic. 3b. Potassium: 


Curve 1. Experimental spectral distribution of photoelectric 
eld in coulomb-calorie, 

Curve 2. Ky X10’, plane (1, 0, 0). 

Curve 3. Ky X10’, plane (1. 1, 0). 

Curve 4. Ky X10’, plane (1, 1, 1). 


body-centered cubic lattices (vp =a*/2) 
km = (62)!/a. 


The lattice constant a=4.24x10-* cm _ for 
sodium and a=5.33X10-* cm for potassium. 
For the work function we take E,=2.29 ev for 
sodium'® and E,=2.01 ev for potassium.!” 
From these data the low frequency limit of light 
absorption and the threshold frequency can be 
calculated from (20) and (23). These and the 
corresponding values of m=g/a are given in 
Table I. The values of », for the non-existent 
g=[(1/a), 0, 0] is also given. The threshold 
frequency depends upon what lattice plane is 
the surface of emission, since the value of g, 
depends upon the direction of x (direction 
normal to the surface of emission). The table 
gives threshold frequency for three different 
lattice planes. For sodium the plane (1, 1, 0) 
gives the lowest threshold 5.9110"; whereas 
for potassium, the plane (1, 1, 1) gives the lowest 
threshold. 5.69 10". These are higher than the 
corresponding ordinary thresholds of emission 
E,/h (5.53X10" for sodium and 4.8410" for 
potassium). The low frequency limit of light 


16 M. W. Mann, Jr. and L. A. DuBridge, Phys. Rev. 51, 


120 (1937). 
17R, Suhrmann and H. Theissing, Zeits. f. Physik 52, 


453 (1928). 
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absorption is 3.14 10" for sodium and 4.96 x 10" 
for potassium. Although experiments do not 
show strong absorption for these metals in the 
visible and ultraviolet regions, but, as pointed 
out by Kronig,'* the “conductivity” (nkyv), 
characterizing the absorption, computed from 
experimental data is 10 times larger than the 
value calculated theoretically without taking 
into account the photoelectric absorption. So, 
although light is weakly absorbed in these 
regions, it should be almost entirely due to such 
absorption. 


APPROXIMATE ESTIMATION OF THE 
PHOTOELECTRIC YIELD 


Without data on electron wave function, the 
photoelectric yield cannot be calculated. We 
shall attempt to estimate approximately its order 
of magnitude. Using (16) we get from (5) 


2 


h 
vH= | 
2rm 





(rg*+2-k) in| 


2 2 


= Vice-k)= 
ae 








g. 


Substituting this into (15) and writing 
(| Ci,e+229|/?) for the average value, we get for 


normal incidence 
=( € ) 1 2xm 
vy? \2mc/ 82° hh? 








4c 
g=1.395x10-| 
T 





xl cal eS fence} 
coulomb 
XK | Ce, e4220 |") ———— 
calorie 
: coulomb 
= KK | Ci, e229 |?) ———— (24) 
calorie 


all factors being in e.s.u. The integrals {/ f'dudv 
can be calculated from Eq. (22). From the data 








TABLE II. 
Sodium 
’ 10" 4.96 6.0 7.0 9.0 11.0 13.0 15.0 
Ry 107" 0 0.666 0.889 1.583 1.597 1.081 0.881 
Potassium 
> 104 3.14 5.2 6.0 7.0 8.0 9.0 11.0 11.82 
Ry 10 0 0.083 0.274 0.417 0.624 0.842 1.80 2.61 
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$1 





of Ives and Briggs’ for m and k, the factor K, 
is calculated and shown in Fig. 3. The value of 
JS JS dudv depends upon what lattice plane is the 
emission surface; K, is calculated for three dif- 
ferent planes, i.e., (1, 0, 0), (1, 1, 0), and (1, 1, 1). 
The rate of consumption of energy per unit 
volume by the excitation of electrons is 


P,=(N/Q)hy. 


The number of excited electrons per unit time 


TABLE III. 











ts 3.21 3.33 3.62 3.88 


Earg—Eo 0.64 0.60 0.50 0.43 
Co, 289 0.2402 0.2214 0.1720 0.1407 

Eseg—Eo 2.27 2.10 1.79 1.54 
Co, 2"9 0.2718 0.2565 0.2317 0.2142 

Earg—Eo 4.84 4.55 3.95 3.47 
Co, 29 0.2045 0.1822 0.1568 0.1513 








N is given by (7) 


2h e Fame 
P,= |» —(-—) = “ff dudo| 
82° \ 2mc 
siittesiaililiath (25) 


The rate of loss of energy per unit volume by the 
light wave is 





2 


2r 
P= }nky|e|?=—nkv*| @|?. (26) 
a 


The ratio is 


Pr C 2h 2 ae 
=-15 Ef fous 
P, 0 & 


gti ioy (27) 


The integrals §fdudv can be calculated by 
using (19). As pointed out above, the energy loss 
of the light wave is almost: entirely due to 
photoelectric absorption, so the ratio P,/P; 
should be nearly equal to one. Table II gives cal- 
culated values of R,. We see that to make 
P./P™\1 the factor (|Ck.e+2e9|*) should be 
~10". 

Although the factor (| Ci,2+2%9|")w is not the 
same in (27) and (24), as it is not averaged over 
the same area in the two cases, its magnitude 








2mc 











%R. de L. Kronig, Nature 133, 211 (1934). 


19H. E. Ives and H. B. Briggs, J. Opt. Soc. Am. 26, 238 
(1936) ; 27, 181 (1937). 
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will not be very different. We have stated 
that / should be ~10-7 cm, so the product 
K | Ci,e+229|?)w~107. From curves K, in Fig. 3, 
we see that the photoelectric yield will be 
between 10-* to 10-* coulomb/calorie. The ex- 
perimental values of absolute photoelectric yield 
for these metals are not very reliable because of 
the difficulty of obtaining perfectly clean surfaces 
for these very active metals. In Fig. 3 the curves 
of Hill? (computed from data of Mann and 
DuBridge'®) for sodium, and Suhrmann and 
Theissing'’ for potassium are reproduced. It is 
seen that our estimated photoelectric yield is 
comparable with the experimental values in 
order of magnitude. It may perhaps be too-small 
for sodium and too large for potassium. Anyway, 
for both metals the photoelectric yield calculated 
from surface effect is the same in order of mag- 
nitude as our estimated values.” As to the shape 
of spectral distribution curves, it cannot be 
predicted in our case without detailed knowledge 
of electron wave functions, since (| Cx, x+229|?)a is 
a function of frequency and the curves K, do 
not, therefore, represent the shapes of the dis- 
tribution curves. Furthermore, ‘the metal is 
usually not a single crystal, and the surface is 
not a single plane; therefore the spectral dis- 
tribution curve should be a combination of those 
for different planes. 
Seitz gives for lithium the values of 


oy, 
Co, on™ J Yiuar 
Ox 


for several values of g. In Hartree’s units 
the results are” those shown in Table III. 
We see that the larger the radius of atomic 
spheres 7,, or rather the larger the part of the 
atomic volume in which y is constant, the smaller 
becomes Co, 279. The radius of atomic spheres for 
sodium is 3.9 (Hartree’s unit), and the region 
around the nucleus, within which the wave 
function is not constant, is about the same as 
for lithium.”* Therefore Co,2,, should be 0.1 to 
0.2 in Hartree’s units. The same is probably 
true also for potassium which has still larger 


” A. G. Hill, Phys. Rev. 53, 184 (Fig. 8), (1938). 

21 R. E. B. Makinson, Proc. Roy. Soc. Al62, 367 (1937); 
R. J. Maurer, Phys. Rev. 57, 653 (1940). 

2F. Seitz, Phys. Rev. 47, 404 (1935), Table III. 

*3 Cf. reference 13, Fig. 1, and reference 22, Fig. 3. 
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radius of atomic spheres, r,=4.93. We have 
estimated 


(| Ci, %429|7)w~10" absolute units 
= 0.0073 Hartree’s units 
((| Cy, e+2eg!2)w)*~0.085 Hartree’s units. 


Although Co,2r, and ((| Ci, +229|?)a)? are not the 
same, we may judge from this that our estimated 
value is quite reasonable. It is satisfying that 
our equation of absorption (25) with a reasonable 
value for (|Ck,+229|")w should agree with the 
absorption given by experimental values of 
and k [(P./P)=1]. 


CONCLUSION 


Although emission due to volume effect does 
begin at a higher frequency, at least for sodium 
and potassium it comes in over most of the 
region for which spectral distribution has been 
investigated experimentally, long before the 
maximum in the distribution curve is reached. 
The estimated order of magnitude of photo- 
electric yield is the same as that calculated from 
surface effect and is comparable with that ob- 
served experimentally. Therefore the volume 
effect cannot justifiably be neglected in the 
theory of photoelectric emission except for the 
immediate neighborhood of the threshold fre- 
quency. It is then worth the labor to calculate 
accurately the spectral distribution curve by use 
of the electron wave functions. It would also be 
interesting to calculate accurately the energy 
loss due to photoelectric absorption to compare 
with the light absorption computed from ex- 
perimental values of m and k. 

Recently Cashman and Bassoe* working with 
barium, which has more stable surface than the 
alkali metals, have observed an abrupt second 
rise in the spectral distribution curve which they 
attribute to the onset of the volume effect. 
Unfortunately, the free electron approximation 
for energy cannot be applied to this metal with 
justification, and we cannot calculate the 
threshold frequency by the same method used 
here for sodium and potassium. We have to wait 
for the solution of electron energy levels for this 
metal. 


% J. Cashman and E. Bassoe, Phys. Rev. 55, 63 (1939). 
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Experiments on Counters with Grids 


SERGE A. Korrr 
New York University, New York, New York 
June 25, 1945 


ONTINUING the investigation of triode counters 
originally undertaken in collaboration with W. E. 
Ramsey,! some further tests have been carried out. In the 
present experiments, several Geiger counters with grids of 
different dimensions were studied. A resistance serving as 
a voltage divider was connected across the high voltage 
source, and the lead from the grid ran to an adjustable 
point on the voltage-divider. In this manner, the field 
between the grid and cathode could be controlled. A Lucite 
window was provided so that a stream of beta particles 
could be projected into the counter. A glass bead, 3 mm 
in radius, was located at the middle of the central wire in 
order to render the device capable of coincidence counting. 
The pulses were examined on an oscilloscope. (See Fig. 1.) 
Voltage pulses corresponding to a simultaneous discharge 
of both halves of the counter were about twice the ampli- 
tude of those corresponding to a discharge of either section 
by itself. An identical counter, with bead but no grid, 
was attached to the same gas-distributing system so that 
the kind and pressure of gas in the two would be identical, 
and the effects produced by the grid could be studied. 

It will be recalled that the effect of the grid is to alter 
the potential distribution inside the counter. The cumu- 
lative ionization (Townsend avalanche) process in the 
Geiger counter discharge takes place in the immediate 
vicinity of the central wire. The rest of the counter serves 
as a volume in which the primary ionization may be 
produced. The customary high counter voltage is necessary 


























Fic."1. Connections for triode counter. 
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to produce the requisite high field near the wire, but the 
electrons can be accelerated toward the high field region 
by a comparatively low field in the outer portions of the 
counter between the grid and the cathode. Thus the grid 
permits a substantial reduction in the over-all voltage. 

One counter was 31 cm long and 7.5 cm in diameter, 
with a 3-mil tungsten central wire and a grid 7.5 mm in 
diameter made of bronze screening. It was filled to a 
pressure of 2 cm of alcohol plus 8 cm argon. It operated 
as a Geiger counter at voltages between 900 and 970. The 
similar counter without a grid and attached to the same 
vacuum system, required 1370 to 1500 volts. Thus, while 
operating voltages are reduced, the width of the plateau is 
also somewhat reduced, but represents approximately the 
same percentage of the total voltage. In an attempt to 
get still lower voltages, a grid 2 mm in radius was tried. 
The counter with this grid, operated at 550 volts when 
filled with a mixture of 1.5 cm of alcohol plus 7 cm of 
argon. Of this 550 volts, 20 volts was the potential differ- 
ence between grid and cylinder. 

Preliminary experiments indicate that the resolving time 
of the counter employing the grid is somewhat less than 
that of the identical one without the grid. The reason for 
this presumably is that the positive ions need to move 
only the small distance to the grid before the field between 
grid and central wire returns to normal. The studies are 
being continued. 


1S. A. Korff and W. E. Ramsey, Rev. Sci. Inst. 11, 267-269 (1940). 





A Theorem on Functional Determinants 


BENJAMIN LIEBOWITZ 
350 Fifth Avenue, New York, New York 
June 25, 1945 


N my paper on “A Calculus of Finite Precision’ | 
gave, without proof, the following theorem: 


—(1/Q)(dQ/ds) = -t. 
Here Q is given by the functional determinant 
ty te see é. 


a 
0a, 0a, Oa; |. 


(1.6) 


Ox OX bea OXn 
Odn-1 OGn-1 O4n~1 








The x; (i=1, 2---m) are functions of n—1 parameters a; 
and of the variable s which measures the distance along a 
curve represented by the x; in the n-dimensional space. 
The vector with components & =dx;/ds is the unit tangent 
to the curve. By differentiation, the a; may be eliminated 
so that the & may be expressed as functions of the x; alone. 

The proof of Eq. (1.6) for three dimensions or less is 
simple and straightforward; for more dimensions, my 
original (unpublished) proof was also straightforward but 
laborious. I have since found a very simple proof which, 
it is felt, should be published for the sake of completeness. 

To prove Eq. (1.6), rewrite the determinant in the form 


1/Q=|a,'|, 
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where 
at =tk ax! = Ox%/Oat-1 (l# 1). 


By the formula for differentiating determinants 


(d/ds)(1/Q) = 2A *(d/ds)ax! 


where A;* is the co-factor of a’. Now 


d dt, Ot, dXm ,0a,! 

Say) me ot om Dt Sm BOE Gt, 

ds° GS mOin GS m OSm 
Also, for 11, 


d 


dq _@ OX, = 9 dx _ Atk AXm _ yAax' , 
ds" dsdai-1 dai. ds 


= a= =a, —t 
mOXm OAi-1 m OXm aie 





provided that the x; are uniformly continuous in s and the 
a;. Hence, for all values of j, 


‘. = dan) i 
oo 
d 1 y Oak 1 0a; 
S38 8 F Apafe an B15 eet 
sO mm om ax km QQ OXm 
1 yay! 1 
Qv ax Q 


1 B. Liebowitz, Phys. Rev. 66, 343 (1944). 





Detection of Metastable Ions with the 
Mass Spectrometer 


J. A. HippLe AND E. U. Conpon 
Westinghouse Research Laboratories, East Pitisburgh, Pennsylvania 
June 4, 1945 


T has usually been customary to operate mass spectrome- 
ters with the filament and electron gun at a positive 
potential with respect to the analyzer section which is 
grounded. Under these conditions the ions pass through the 
exit slit and strike the ion-collector plate with the full 
energy acquired in the ion gun since this plate is at ground 
potential. Recently with an instrument of the sector type 
employing a magnetic deflection of 90°, the first slit in the 
ion gun was connected to ground while the analyzer section 
was made negative with respect to ground. Since the ions 
are formed in a region one or two volts above ground, they 
reach the collector plate in this arrangement with an energy 
of only a few electron volts. 

This change has had a very interesting effect on the 
spectra of hydrocarbons. In Fig. 1 a portion of the spectrum 
of normal butane is shown. The lower curve shows the 
baseline with the conventional arrangement, and the upper 
curve shows the improved baseline when the ion source is 
grounded. These curves were obtained with an automatic 
recorder employing a non-linear scale with provision for 
recording large peaks (such as mass 43 which was allowed 
to go off scale since it is of no concern at the moment). The 
“hump” on the side of the peak at mass 39 has not previ- 
ously Deen investigated carefully but was usually attributed 
to the formation of C;H;* with kinetic energy. However, it 
should still be present in the upper curve of Fig. 1 if this 
explanation is correct. Its disappearance means it must be 
associated with a loss in energy. Similarly, the disappear- 
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ance of the diffuse peaks around 30.5 and 32 requires a loss 
in energy. 

The experimental results may be explained by supposing 
that some of the ions dissociate after acquiring energy in the 
ion gun, i.e., if metastable ions are formed in the ion source, 
If this dissociation occurs in the region of the instrument 
between the ion source and the magnet, the ion will have 
an apparent mass m* which has the following relation to 
mo, the mass of metastable ion, and m,, the mass of the ion 
to which the metastable dissociates, 


m* = m,?/mo. 


This simple relation assumes that the metastable ion 
dissociates into two fragments with negligible release of 
internal kinetic energy. In general, of course, there may be 
some energy release. Also the dissociation processes in 
general will follow an exponential decay distribution in 
time, and so some will occur in the magnetic analyzer and 
other parts of the instrument giving rise to the observed 
unsharpness of the peaks. On dissociation the ion has 
therefore only that fraction of the kinetic energy repre- 
sented by its fraction of the original mass and is therefore 
unable to reach the collector when there is a full retarding 
field between the main analyzer tube and the collector. 

The peak around mass 32 can be explained by the 
dissociation 

C,Hi9t—>C;H7*+CHs, 
giving 
m* = (43)?/58 = 31.9. 
In this case the parent ion dissociates, and one fragment is 
the most abundant ion in the n-butane spectrum. 

Similarly the ionization around mass 30.5 in the spectrum 

can be explained by the dissociation 


C,H 10° —>C;Het +CHy, 
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Fic. 1. Portion of the spectrum of normal butane obtained on an 
automatic recorder employing a non-linear scale. The two curves show 
the effect on the spectrum of operation of the mass spectrometer with 
and without the ion source grounded. The disappearance of the diffuse 

k at mass 32, for instance, is attributed to the presence of metastable 
ions in the mass spectrometer. The small residual peak at mass 32 is 
caused by O+impurity in the tube. The original charts have been 
retouched for reproduction. 
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giving 
m* = (42)2/58 = 30.4. 


The hump on the side of mass 39 might be explained by 
m* = (40)?/41 = 39, 


m* = (41)?/43 = 39.2. 


Similar peaks have been noticed in the spectra of a great 
many hydrocarbons, and in al! cases the explanation ap- 
pears to fit into the pattern outlined above. Some other 
observations on mass spectra of hydrocarbons, such as 
variation of relative peak heights from one instrument to 
another, and with initial temperature of the gas, find plausi- 
ble interpretation as effects due to dissociation of meta- 
stable ions. 

We wish to thank H. A. Thomas and R. E. Fox for their 
assistance in this experiment. A more complete report on 
metastable ions in the mass spectrometer is in preparation. 





Further Work on the Artificial 
Production of Mesotrons 


GERHART GROETZINGER AND LLOYD SMITH 
University of Illinois, Urbana, Illinois 
July 3, 1945 


S reported in an earlier publication, one of us observed 

a new penetrating radiation in the neighborhood of a 
cyclotron a few years ago.' A detailed investigation showed 
that this neutral radiation definitely had properties differ- 
ent from those of neutrons and gamma rays.'? The 
suggestion was made that the radiation consisted of low 
mass uncharged mesotrons, a hypothesis which was 
strengthened by the reported occurrence of charged meso- 
trons of from 10 to 20 electron masses in cosmic radiation.* 
Recently, cloud chamber pictures have been reported‘ 
which are claimed to confirm the existence of low mass 
neutral mesotrons. These pictures show electron pairs of 
high energy emerging in opposite directions from an 
aluminum plate placed in a cloud chamber situated some 
distance from the cyclotron and in line with the beam 
striking the target. 

In our earlier publication,' it was stated that the elec- 
trons which caused the coincidences were of energy greater 
than 5 Mev. With the arrangement described at that 
time, experiments have been made to investigate the 
properties of these electrons by interposing slabs of alumi- 
num between the two coincidence counters. Figure 1 shows 
the dependence of the coincidence counting rate upon 
thickness of aluminum (including the equivalent thickness 
of the two counter walls), when a Mn target is bombarded 
by the cyclotron. It is seen that the counting rate decreases 
gradually at first and then quickly drops to zero at a 
thickness of about 12 mm. 

If this behavior is interpreted in terms of a process 
which leads to the production of single electrons only, then 
these electrons must be nearly all of the same energy. The 
second scale of abscissae represents the energy of an elec- 
tron which is just able to penetrate the corresponding 
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thickness of aluminum, showing that the energy of this 
homogeneous group would be approximately 6.5 Mev. 

If, on the other hand, the much more probable process 
of the decay of the neutral mesotron into a pair of electrons 
is assumed, some of the coincidences are caused by one of 
the pair electrons penetrating both counters and the rest 
by pairs created between the counters, one electron tripping 
each counter. From a certain thickness of aluminum on, 
all the counts will be due to mesotrons decaying between 
the counters. The third energy scale gives the energy of 
one of a pair of electrons which could just penetrate half 
of the corresponding thickness of aluminum having been 
created in the center of the slab. If this second interpreta- 
tion is correct, the rest energy of the mesotrons is 2(3.5+.5) 
=8 Mev=16 me’. 

An estimate of the cross section for the production of 
neutral mesotrons is now possible, since the energy of the 
electrons determines the thickness of lead near the counters 
which is effective as an electron producing layer and thus 
determines the effective solid angle subtended by the 
counter arrangement at the target. If we assume that all 
counts recorded behind 19 cm of lead! are due to the decay 
of neutral mesotrons which have been slowed down or 
stopped in the lead, the counting rate behind 28.5 cm of 
lead gives a measure of the rate of absorption of the 
mesotrons and enables one to calculate the number of 
mesotrons leaving the target. Because of the small mass of 
the particle, the forward direction is probably preferred in 
the emission process; if we estimate that mesotrons leave 
the target in the forward direction ten times more often 
than the number per unit solid angle averaged over all 
directions, the total cross section for production is found 
to be 10-** cm*. The mechanism by which the particles 
are created is probably similar to the “bremsstrahlung” 
process postulated in the case of cosmic ray mesotrons, 

We have also performed an experiment on the production 
of charged mesotrons by 11-Mev deuterons using a method 
of delayed coincidences. For this work, a portion of the 
deuteron beam was deflected in a magnetic analyzer and 
piped through the water tanks which surround the cyclo- 
tron.’ The beam was made to strike foils of Ail, Ni, Cu, 
and Pb just thick enough to stop the deuterons. In line 
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Fic. 1. Coincidences per 10"’ deuterons striking the target as a function 
of thickness of aluminum between the counters. 











56 LETTERS TO 


with the incident beam and separated from the target by 
a 15-cm long channel formed by lead blocks was placed a 
counter telescope consisting of two thin-walled counters. 
Behind the telescope was an aluminum plate inclined at 
an angle of 45 degrees to the direction of the telescope. 
Beside this plate, out of line with the target and telescope, 
was placed a third thin-walled counter. The purpose of the 
arrangement was to detect charged mesotrons which were 
created in the target and stopped in the aluminum plate 
after traversing the telescope, the decay electrons tripping 
the third counter. The circuits were adjusted so that 
processes were recorded in which the third counter was set 
off 1.0 to 1.4 microseconds later than the two counters of 
the telescope. With a beam of 10~* ampere at the target, 
real delayed coincidences occurred at a rate of two per 
minute. 

However, the counting rate was not changed appreciably 
by removing the aluminum plate or by interposing several 
millimeters of lead between the telescope and the delay 
counter to shield the latter completely from the expected 
decay electrons. Therefore, the observed counts can not 
be attributed to charged mesotrons decaying into electrons. 

If charged mesotrons of low mass are actually created 
when the deuterons strike the target, the failure of this 
experiment may be laid to a combination of three causes: 

(1) Production cross section too small. If the lifetime of 
these particles is of the order of 10~* second, a cross section 
for production greater than 10-* cm? would be sufficient 
for their detection in this experiment. 

(2) Lifetime too long or too short. If the cross section 
for production is 10-** cm?, as for the uncharged mesotrons, 
then particles of a lifetime longer than 10-5 second or 
shorter than 5X 10~7 second would pass unnoticed. 

(3) Kinetic energy acquired by the mesotron at the 
moment of creation too small for it to penetrate the target 
material and the counter walls. A charged particle of rest 
energy 20 mc? and a kinetic energy of 1 Mev or less would 
easily be stopped before reaching the aluminum plate. 


1G. K. Groetzinger, P. G. Kruger, and L. Smith, Phys. Rev. 67, 52 


ge). 
2G. K. Groetzinger and L. Smith, Phys. Rev. 67, 53 (1945). 
3M. Schein (private eqaegeneten. 
4P. G. Kruger and L. W. Smith, abstract presented before the 
Stanford meeting of = American Physical Society, July 14, 1945. 
5P. G. Kruger, G. K. Groetzinger, ef al., Rev. Sci. Inst. 15, 333 


(1944). 





Refraction Effects in Electron Diffraction 


LORENZO STURKEY AND Lupo K. FREVEL 
The Dow Chemical Company, Midland, Michigan 
July 3, 1945 


iy an attempt to use MgO smoke as a calibration 
specimen for electron diffraction patterns in a camera 
with a 72.3-cm specimen-to-plate distance, rather unex- 
pected diffraction patterns were obtained: 

(a) The (hhh) reflections appeared as resolved doublets. 

(b) The (h00) reflections were very sharp. 

(c) Other reflections were considerably broadened. 
(This phenomenon for MgO has also been observed by 
J. Hillier and reported without explanation in’a contributed 


THE EDITOR 


paper at the American Physical Society meeting at Ohio 
State, June, 1945. Hillier used a camera with a 20-cm 
specimen-to-plate distance.) 

Since x-ray back reflection patterns taken with CuKa- 
radiation of MgO smoke deposited on a glass microscope 
slide showed no signs of deviation from the usual cubic 
structure, it must be concluded that the broadening and 
doubling of reflections other than (00) are characteristic 
of electron diffraction. It is well known that electron 
diffraction patterns are very sensitive to the shape of the 
diffracting particles. The observed splitting of the non- 
pinacoid reflections is attributable to the uniform cubic 
habit of MgO smoke. Electron rays reflected from the 
(h00) planes pass through the crystal as through a parallel 
plate and suffer no appreciable change in direction other 
than the change by diffraction; but rays reflected from 
prism planes often pass through an effective prism and 
suffer an additional change in direction due to refraction 
at the crystal faces. To check this deduction, the following 
observations were made: 

(1) CdO smoke, which also shows a pronounced cubic 
habit, gives electron diffraction patterns exhibiting the 
same peculiarities as MgO smoke. 

(2) MgO was prepared in irregularly shaped particles 
ranging in size from 200A to 500A as observed in the 
electron microscope. Electron diffraction patterns of this 
material showed no anomalies, and all diffraction rings are 
essentially equally sharp, though not so sharp as for MgO 
smoke. 

(3) Diffraction patterns of MgO smoke were taken at 
various electron accelerating voltages. It was frequently 
observed that, at lower voltages, the center of the broad- 
ened (220) reflection decreased in intensity. Since this 
central part is due to diffraction through a consistently 
thicker portion of the crystal than the side parts, absorp- 
tion tends to decrease the intensity at lower voltages. 

(4) ZnO smoke patterns show broadening of the (100) 
and (101) reflections. This material has a spike-like habit 
instead of cubic. 

For a cube-shaped crystal, the deviation of the electron 
rays in the diffracted direction—the only direction affecting 
ring broadening—may be approximately calculated by 
finding an “effective refracting triangle’ for the ray. This 
effective triangle is found by passing a plane through the 
cube perpendicular to the diffracting plane and including 
the incident and reflected rays from this diffracting plane. 
The intersection of this plane with the diffracting plane is 
the base of the effective triangle, and the intersections with 
the faces of the cube on the “reflecting” side of the dif- 
fracting plane form the other two sides. The deviation 6 
in radians from the Bragg angle is given approximately by 
the following expression: 
6=(+[tan* r,+sIin 7; cos r:] 

+[tan? re+sin r2 cos re }(P/2E), 
where 
r, and rz are the angles of refraction at sides of triangle, 
sign of brackets being megative if face normal is 
between ray and diffracting plane, and positive 
otherwise 
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P=inner potential of crystal 
E=accelerating potential of electrons. 


A calculation of the inner potential of MgO from the 
doublet spacing of the (111) reflection gives a value of 
12+4 volts, which agrees in order of magnitude with 
measurements, etc., of inner potentials for other ionic 
crystals." 

A more detailed analysis of the effect of refraction on 
electron diffraction patterns will be published later. 


1(See Thomson and Cochrane, Theory and Practice of Electron 
Diffraction, Chap. X.) 





K-Capture and e+-Emission of 6.7-hr. Cd: An 
Experimental Verification of the 
Fermi Theory of $-Decay 


H. Brapt, P. C. GuGetot, O. Huser, H. Mepicus, 
P. PREISWERK, AND P, SCHERRER 


Physikalisches Institut der Eidgendssischen Technischen 
Hochschule, Zurich 


July 3, 1945 


T may be possible to discriminate between different 

theories of 8-decay by the determination of the relative 
probabilities of K-capture and e*-emission. The experi- 
mental difficulties of an exact determination of this kind 
are evident; until now, no exact value of the quotient 
Q=px/p+ for an allowed transition was known. 

We found, that the 6.7-hr. Cd!®7/10® which is mostly 
decaying by K-capture, emits also positive electrons. The 
e*-spectrum was measured with the magnetic spectrometer 
and the upper limit of the spectrum found to be 
E=(320+10) kev. The corresponding annihilation-radia- 
tion was detected by (y-y)-coincidences, following the 
6.7-hr. period. 

Both disintegrations, K-capture and e*+-emission, lead to 
a metastable Ag'®’/!*-nucleus which goes over to the stable 
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Ag with a half-period of (44.3+0.2) sec. Only 1 percent of 
this transition is connected with the emission of a 93.5- 
kev y-ray; in 99 percent, according to the high internal 
conversion coefficient, a K, L, or M electron is ejected. 
Therefore, the number of conversion-electrons of the 93- 
kev radiation, emitted by the 44-sec. Ag!®7/! in equilibrium 
with the 6.7-hr. Cd, is practically equal to the total number 
of the Cd-disintegration-processes. The relative proba- 
bilities of K-capture and e*-emission are therefore in this 
special case easily determined, being the proportion of 
the total number of conversion-electrons to the total 
number of positrons from the same source of 6.7-hr. Cd. 

The experimental value, determined with the magnetic 
spectrometer is found to be 


Qexp = Px/Ps =320+20 


in excellent agreement with the theoretical value calculated 
according to the Fermi theory of 8-decay. 


Qtheor Fermi = 342. 


The Konopinski-Uhlenbeck modification leads to a value 
60 times greater than the one actually found: 


Qxu = 19800. 


With A,=8.97X10-* sec! and W.,/mc?=1.62 the 
matrix element of the Fermi theory turns out to be 
M=0.125, the e*-decay therefore corresponding to an 
allowed transition. 

Besides the radiation of 0.51 Mev resulting from annihi- 
lation of the e+, we were able to measure the energy and 
intensity of the weak nuclear y-radiation of the 6.7-hr. Cd. 
The energy of this y-radiation is 0.842 Mev and the abso- 
lute intensity is 0.042 quanta per disintegration. The 
K-capture leading to the excited level of Ag with 0.842 
Mev energy is, therefore, at once forbidden. A detailed 
report will be published in the Helvetica Physica Acta, 
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A meeting of the Metropolitan Section of the American Physical Society was 
held on Saturday, June 2, 1945, at the Polytechnic Institute of Brooklyn, 85 
Livingston Street, Brooklyn, New York, as a joint local symposium with the 
Polytechnic Institute and the Society of Rheology on the Ultracentrifuge in 
Polymer Research. The program was as follows: 













10:00 a.m. 
Introductory Remarks, J. B. NicHoits, Chairman of the Symposium (E. I. duPont 
deNemours & Company, Wilmington, Delaware). 











10:30 a.m. 
Special Problems in Ultracentrifuge Research on High Polymers, C. O. BECKMANN 
(Columbia University, New York City). 


1:00 p.m. 
Business Meeting of the Metropolitan Section of the American Physical Society. 













1:30 p.m. 

The Size Distribution of Gelatin, GEorGE SCATCHARD (Massachusetts Institute of 
Technology), J. L. ONctey (Harvard University), AND J. W. WiLLiAMs (University 
of Wisconsin). 


3:00 p.m. 
Studies on the Molecular Weight and Shape of Linear Polymers, Paut Dory, 
SEYMOUR SINGER AND Kurt G. STERN (Highpolymer Research Bureau, Polytechnic 

Institute of Brooklyn). 


4:00 p.m. 
General Discussion. 


























At the business meeting of the Metropolitan Section the following officers 
were elected for the season 1945-1946: 











Chairman, H. Mark 
Brooklyn Polytechnic Institute 
Vice Chairman, H. A. Barton 
American Institute of Physics 
Secretary- Treasurer, W. H. Crew 
National Defense Research Committee 









Elected Members of the Executive Committee: 









Edith Quimby, Dept. of Radiology, Columbia University 
J. B. Fisk, Bell Telephone Laboratories 








A vote of thanks was tendered to W. S. Gorton for his service of twelve years 
as Secretary-Treasurer of the Metropolitan Section. 
W. S. Gorton, Secretary-Treasurer 
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